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ABSTRACT 


The  Navj’  Remote  Ocean  Sensing  System  (  N-ROSS  )  satellite  is  being  developed 
to  supply  accurate  data  on  ocean  parameters  for  fleet  operations.  A  Low  Frequency 
Microwave  Radiometer  (  LFMR ),  a  large  flexible  reflector  attached  to  an  angled 
flexible  boom,  is  a  sea  surface  temperature  sensor  on  this  satellite  which  rotates  at  15 
RPM.  The  dynamic  interaction  betw'ecn  the  reflector  and  the  boom,  and  the  effects  of 
the  reflector  orientation  and  flexibility  on  the  pointing  error  of  the  LFMR  during  a 
spin-up  procedure  are  investigated  by  performing  dynamic  simulations.  Dynarrucal 
equations  of  this  flexible  multibody  system  are  formulated  using  Kane's  method. 
Efficient  computer  simulations  were  achieved  by  developing  a  FORTRAN  program 
and  using  Dynamic  Simulation  Language  (  DSL  ). 
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A.  BACKGROUND 

The  Navy  has  great  interest  in  obtaining  accurate,  real-time  data  on  various  ocean 
parameters.  Presently,  research  is  being  conducted  on  the  merits  of  using  a  satellite  system 
to  collect  this  data.  The  Navy  Remote  Ocean  Sensing  System  (N-ROSS)  satellite  has  been 
developed  to  achieve  this  goal  [Ref.  1].  This  satellite  scans  the  earth’s  surface  and  will 
provide  the  Navy  with  such  parameters  as  wind  speed,  wind  direction,  ocean  temperature, 
ice  edge  detection,  wave  height,  and  ocean  photography  [Ref.  2]. 

One  of  the  sensor  systems  to  accomplish  this  goal  is  the  Low  Frequency  Microwave 
Radiometer  (LFMR).  The  LFMR  consists  of  a  large  reflector  dish  connected  to  an  angled 
boom.  The  LFMR  scans  the  ocean  surfaces  and  records  the  temperature.  Increased  scan¬ 
ning  area  is  accomplished  by  spinning  the  reflector  at  15  RPM.  A  rigid  electronics  box  is 
attached  at  the  base  of  the  LFMR  boom.  Therefore,  the  rotating  LFMR  system  consists  of 
the  flexible  boom  and  reflector  and  the  rigid  electronics  box.  Considering  the  nature  of  the 
system,  accurate  analysis  and  suppression  of  the  vibration  and  deflection  of  the  LFMR  is 
essential  for  the  satisfaction  of  the  pointing  error  requirement  of  the  system  and  for  the 
accomplishment  of  the  mission  of  the  N-ROSS  satellite[Ref.  2]. 

B .  PROBLEM  STATEMENT 

A  three-dimensional  dynamic  analysis  has  been  performed  using  a  model  that  consists 
of  a  flexible  angled  boom  with  a  point  mass  at  the  tip  [Ref.  3].  Lagrangian  equations  of 
motion  were  employed  to  achieve  a  good  dynamic  simulation. 


In  this  study,  the  LFMR  is  analyzed  as  a  multibody  system.  A  multibody  approach 
will  provide  a  means  to  analyze  the  dynamic  interaction  between  substructures  of  a  compli¬ 
cated  structure  with  greater  ease.  In  dynamic  analysis  of  a  multibody  system,  maintaining 
computational  efficiency  is  a  key  factor  and  “The  emphasis  of  researchers  working  with 
multibody  systems  has  been  the  expanded  generality  of  mathematical  models  and  the  for¬ 
mulation  of  equations  of  motion  that  are  amenable  to  computer  solution.”  [Ref.  4] 

Therefore,  it  is  the  intent  of  this  thesis  to  apply  an  efficient  multibody  dynamic  analy¬ 
sis  to  the  LFMR  system.  The  LFMR  is  broken  down  into  two  bodies;  the  first  body  is  the 
flexible  boom  and  the  second  body  is  the  flexible  reflector.  Kane’s  method  is  used  to  for¬ 
mulate  these  equations  of  motion,  “since  it  combines  the  computational  advantages  of  both 
Newton’s  laws  and  the  Lagrangian  formulation.”  [Ref.  4] 

C .  THESIS  OUTLINE 

In  Chapter  II,  the  analytic  model  of  the  LFMR  used  for  this  thesis  is  described.  The 
dynamical  equations  are  formulated  using  Kane’s  method.  This  method  is  briefly 
explained  prior  to  the  formulation. 

Chapter  III  contains  an  explanation  of  the  three  computer  programs  used  to  solve  the 
equations  of  motion.  NASA  Structural  Analysis  (NASTRAN)  calculates  the  mode  shapes 
and  natural  frequencies  for  both  bodies.  A  FORTRAN  program  reads  the  data  output  from 
NASTRAN  and  calculates  the  time  constants.  The  time  constants  are  applied  in  a  Dynamic 
Simulation  Language  (DSL)  program  to  solve  the  simultaneous  differential  equations. 

In  Chapter  FV,  the  simulation  results  are  presented  to  investigate  deflections  and  rota¬ 
tions  at  various  points  on  the  LFMR.  Three  parameters  arc  varied:  (1)  the  flexibility  of  the 
reflector;  (2)  the  in-plane  orientation  of  the  reflector,  and  (3)  the  out-of-plane  orientation  of 
the  reflector.  The  effects  of  these  changes  are  compared. 
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Chapter  V  presents  conclusions  and  recommendations  for  application  and  extension 


of  this  work. 


IL  FQRMIILATIDN  QF-EQUATIQISS-QF  MOTION 

A.  MODEL  DESCRIPTION 

The  configuration  of  the  LFMR  is  shown  in  Figure  2.1.  The  LFMR  consists  of  a 
boom  mounted  on  a  rigid  electronics  box.  A  complex  refiector  is  attached  at  the  tip  of  the 
boom.  The  boom  hinge  is  to  be  rigid  after  the  deployment  of  the  LFMR,  allowing  no  rela¬ 
tive  motion  between  the  upper  and  lower  booms.  The  LFMR  system  rotates  at  15  RPM 
about  a  spin  axis  fixed  in  a  local  reference  frame  of  the  N-ROSS  satellite.  The  spin  axis 
always  points  at  the  earth’s  center  while  the  satellite  is  in  its  orbit.  Therefore,  the  gravita¬ 
tional  force  and  the  centrifugal  force  are  assumed  in  equilibrium.  This  observation  allows 
the  spacecraft’s  local  reference  frame  to  be  considered  the  inertial  (Newtonian)  reference 
frame  and  the  LFMR  to  be  in  a  zero-gravity  environment 

The  model  of  the  LFMR  used  in  this  thesis  is  shown  in  Figure  2.2.  The  angled  boom 
was  analyzed  as  one  body  and  the  reflector  as  the  second  body.  Body  A  represents  the 
electronics  box,  body  B  is  the  angled  boom,  and  body  C  is  the  reflector  dish. 

The  31,32,33  coordinate  system  rotates  in  the  inertial  reference  frame,  with  33 
coinciding  with  the  spin  axis.  The  C 1  ,€2.03  coordinate  system  is  the  local  reference  frame 
of  the  reflector  and  its  origin  is  at  the  connection  point  between  the  deflected  boom  and 
reflector  (point  h). 

B .  EQUATION  FORMULATION 

1.  Kane’s  Dynamical  Equations 

The  equations  of  motion  for  the  LFMR  were  formulated  using  Kane’s  method 
[Ref.  5].  This  method  is  a  Lagrange  form  of  D’Alambert’s  principle,  which  is  equivalent 
to  Newton’s  law  cast  into  a  different  form.  Newton’s  law  for  a  differential  element  is 
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(Eqn.  2.1) 
where 


df  -  ai^dm  =  0 


df  -  differential  force 

3*^  -  acceleration  of  differential  element  on  body  K 
dm  -  differential  mass 

Kane’s  method  introduces  the  generalized  active  and  inertia  forces.  The 
dynamical  equations  are 

(Eqn.  2.2)  +  fJ!  =  0  r  =  1 . total  degree  of  freedom 

F^  represents  the  generalized  active  forces  and  F“  is  the  generalized  inertia 
forces.  Ff  is  the  sum  of  the  dot  product  between  the  partial  velocity  and  the  differential 
force  as  shown  in  equation  2.3. 

(Eqn.  2.3)  F^  =  J  Yr  ‘ 

F“  is  the  sum  of  the  dot  product  between  the  partial  velocity  and  the  differential 
inertia  force  as  shown  in  equation  2.4. 

(Eqn.  2.4)  J  '  (-a*^)  dm 

The  partial  velocity  comes  from  the  definition  of  velocity  where 
V 

(Eqn.  2.5)  y  =  X  ^  '^r  . degrees  of  freedom 

r=l 

is  the  time  derivative  of  the  generalized  coordinates. 

For  this  model,  there  are  three  types  of  dynamical  equations.  Equation  2.6  is 
derived  from  the  large  rotational  motion  of  the  system.  Equation  2.7  is  derived  from  the 
small  boom  motion.  The  number  of  equations  of  this  type  depends  on  the  number  of 
modes  (v  ^ )  used  to  characterize  the  boom  motion.  Equation  2.8  is  derived  from  the  small 
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reflector  motion.  The  number  of  equations  of  this  type  depends  on  the  number  of  modes 
(V2)  used  to  characterize  the  reflects  motion. 

(Eqn.  2.6)  F 1  +  F"  =  0 

(Eqn.  2.7)  F^  +  j  +  F,*j  =  0  j=l . Vj 

(Eqn.  2.8)  Fi+v,+k  "*■  F)*V|*k  “  0  k=l,...,V2 

2.  Generalized  Coordinates  and  Generalized  Sneed 

The  first  generalized  speed  u  1  is  defined  as  follows: 

N  A 

(Eqn.  2.9)  U  1  =  iii.  •  83 

N  A 

where  is  the  angular  velocity  of  body  A  in  the  inertial  reference  frame.  Subsequent 
generalized  speeds  are  defined  from  the  generalized  coordinates  C]  j,  cfj  which  are  model 

coordinates  of  bodies  B  and  C.  respectively,  as  follows: 

(Eqn.  2.10)  Ui*i=qi  (i  =  1 . V,) 

Ut+V|  + j  ■  (j  =  1 . V2) 

where  V  ]  and  V2  represent  the  number  of  modes  used  to  describe  the  displacement  of  body 
B  and  C,  respectively.  Thus,  there  arc  1  +  V 1  +  V2  generalized  speeds. 

3.  Aneular  Velocities 

The  angular  velocity  of  body  C  with  respect  to  body  A  is  the  time  derivative  of 
the  small  rotation  of  the  hinge  point  H.  Point  H  is  the  connection  between  body  B  and 
body  C.  This  angular  velocity  can  be  represented  using  matrices  as  follows: 

r  1  '^'3-92’]  f  _0  -'|'3'|'2l  r  0  -'|'3'i'2 

(Eqn.  2.11)  "'(x)  =  [cllc]  =  -'^'3  ’  'I'l  4*3  0  =  M'3  0  -'I', 

"  Lt2"4'i  1  J  L-'I'z'i'i  0  J  L-92  4'i  0  _ 

Vi  .  . 

where  4*  j  =  S  4^ k  (H )  U 1  + and  9  k  the  small  angle  displacement  for  the  i  ^ 
k=l 

degree  of  freedom  and  the  k*  mode  at  point  H.  Therefore,  the  angular  velocity  is 
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(Eqn.  2.12) 


SLi  ♦  'l'2  'l'3  £3 


4. 


Velocity  and  Ao 


£ 


a. 


(Eqn.  2.13) 


Velocity  and  Acceleration  of  Body  A 

From  the  generalized  speed  definidcxi,  die  angular  velocity  of  Body  A  is 
NA 

^  =  Ut  ^ 


The  angular  acceleration  is  merely  the  time  derivative  of  the  angular 
velocity,  because  the  83  axis  is  fixed  in  the  inertial  reference  frame. 


(Eqn.  2.14)  Ci  -  3.3 


The  velocity  of  A  is  defined  relative  to  point  Q,  the  origin  of  the  a  ^  ,02,33 


coordinate  system.  Point  Q  is  Exed  in  the  inertial  reference  frame.  Since  body  A  is  a  rigid 

body,  the  velocity  is  merely  the  rotational  velocity  of  the  mass  center  A*. 

NA**  ,^NA 

(Eqn.  2.15)  V  =  Ciil  =  ”bu,32 

where  b  is  the  distance  from  point  Q  to  the  mass  center  A*. 


Similarly,  the  acceleration  of  body  A  is  the  rotational  acceleration  com¬ 
prised  of  tangential  and  normal  components. 

NA**  NA  NA  NA 

(Eqn.  2.16)  3  =  iii.  (ill  "ba))  +  a  x  ('bai) 

2 

=  bu ,  3i  -  bu.  32 


The  a )  .02.33  coordinate  system  was  chosen  so  that  the  mass  center  of 
body  A  moves  in  plane  motitm,  simplifying  the  equations, 
b.  Velocity  and  Acceleration  of  Body  B 

The  position  of  an  arbitrary  point  P  on  body  B  is  described  in  the 
a  1 ,32,03  coordinate  system. 
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(Eqn.  2.17) 


If  =  Pi  ai 


*  Pz  *  P3  2.3 


S't'i  <1|  ai 

i=l 


'^1  n  7 

S<t>i  ^123 
i=l  i=1 


B  P 

The  position  vectt}r  B  is  the  sum  of  the  undeformed  position  vector  and 
the  summations  of  the  translational  mode  shape  deformations.  P 1 ,  P2>  P3  describe 


the  undeformed  position. 

The  position  of  point  P  changes  with  time,  therefore  the  velocity  of  point 


P  is  the  time  derivative  of  the  position  added  to  the  rotational  velocity. 


(Eqn.  2.18) 


NP  B.P  NA  BP 
V  =  B  iil  ^  S 


Vi  ^  V,  Vi  ^ 

=  Sf  j  Uhl  ai  +  S<|>i  Ui+i  22  +  'Lh  ui+i  23 


i- 

=  1 

i=l 

+  u, 

V 

22  -  u, 

M 

Pi*  E'liqi 

i=1 

i=I 

The  acceleration  of  point  P  is  derived  from  differentiating  the  velocity 


with  respect  to  time. 
(Eqn.  2.19) 


NP  NA 

Ci  =  fii  X 


BP  NA  ,NA  BP.  ^NA  BF 
a  *  iii.  (jlii  ^  B  )  •*•  2iii.  X  £ 


BP  BP 
+  2 


'  -u’l 

Pz* 

2 

-  u, 

Pi*  X4'iqi 

< 

i=l 

i=i 
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•  • 

> 

2u, 

u 

V 

♦  E<l>i  UUi 

i=l 

i=1 

> 

r 

"  ut 

- 1 

y 

+ 

—  X 

Z _ 1 

2 

-  U, 

Pz- 

1 

II 

_ 1 

— 1 

II 

_ 1 

2u 


Si2 


*  S'l>i 
1=1 


c.  Velocity  and  Acceleradonof  Body  C 

The  position  vector  of  an  arbitrary  point  X  on  body  C  is  described  in  the 
local  coordinate  system  C  analogous  to  the  representation  of  a  pmnt  on  body  B. 

/  \  r  ^ 

Vz 

«  II 


(Eqn.  2.20) 


= 


V 


X2*  Ztjqj 
j=i 


V 


y 


X3* 

j=i 


£3 


V  y 

The  velocity  of  point  X  in  the  inertial  reference  frame  is  the  sum  of  the 
velocity  at  the  hinge  point  H  in  the  inertial  reference  frame  and  the  rotational  velocity  with 
the  time  derivative  of  the  position  calculated  in  the  local  coordinate  system  C. 
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(Eqn.  2.21) 
where 


V  =  T .  (S,"  X  -ff )  .  ^ 


=  U,  23  +  'I'l  £l  'i'2  ^  '^'3  ^ 


'a=*  ^t^Wui.rut  H2*2't>i(H)qi  'ai 

i=i  i=l 


I<t>?(H)Uj.,*u,  H,*  X<(>j(H)q 

i=l  i=1 


i=l 


Thus  the  velocity  of  point  X  is 


NX 

(Eqn.  2.22)  y  = 


i=l  i=l 


i=l  «=1 


r  f 

* 'f'2 

i=l  j=l  J-' 

V. 
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/ 

X2*  f 

j=l 


V 

/ 


/ 


j=l 


^2  _ 


Xl+  X+j^j 

j=1 


V 


-'/'t 


y 


X3*  i 
j=l 


/ 


V 


j=l 


X2+  I  jqj 
j=1 


y 

h'^'2 


r£2 


J 

\  r 


\ 


y 


Xi+  Xljqj 
j=l 


r 

\ 

r 

(  M 

V?  ... 

^  -u 

<  u, 

Xl  + 

M 

-e- 
—  > 
jCkI 

►  il,  +  -< 

Ul 

J=1 

)=1 

k. 

) 

k. 

L  p 

r 

? 

\ 

Yz 1 

<  Ui 

X3^ 

kl  •— 

-e- 

^  d.3* 

j=i 

k. 

V 

) 

4 

y 

lil2  ^ 


A  third  nonoithogonal  coordinate  system  (d)  is  defined  for  computing  the  velocity 
and  acceleration  of  body  C  more  easily. 

33  £i  =  SLi 
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The  following  tables  depict  the  relati(»ship  between  the  three  coordinate  systems. 


Ai _ iz _ 4 


ill 

1 

0 

d2 

-1 

-^3 

0 

^2 

0 

(not  orthogonal) 


The  acceleration  of  point  X  with  respect  to  the  inertial  reference  frame  is 
expressed  as  the  hinge  point  H  acceleration  and  the  acceleration  relative  to  the  hinge  point 
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NX  ,  NH  ,  NC  ^  |X  ^  ^  ^  (NX  X  ^  2  ^  ^  « 


(Eqn.2.23)  a  =  fl  ^  fii 
where 


,AC. 


NC  NA  ^(Oi)  )  NA  AC 

a  =1  *  “• 


=  U,  a3  ^  4*1  ill  *  *  '^’3  Sl3  * 

=  Ui  4^2  ila  *^3  il3 


and 


NH  ^ 

a.  =< 


-Ui 


Hz*  S't''|(H)qi 
i=i 


2 

-Uj 


Hi+  2 
i=1 


-2UiI<t>^(H)Uur  I<t>>)^i*l 

i=1 


1= 


a\ 


"  Ui 

Hi*  Xti(H)qi 
i=l 

2 

'U, 

Hz*  X't’i(N)qi 
i=t 

- 

* 

Vl 


Vl 


-2u,I<t>>)ui*rX<t>i(H)uui 
i=l  i=^ 


a2  ^ 


'^1 

S  a3 

i=l 
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The  acceleration  equation  for  body  C  was  too  unwieldy  for  higher  terms 


to  be  considered.  Therefore,  only  first-order  terms  were  used  to  produce 


(Eqn.  2.24) 


NX 

a  = 


" 

(  \ 

NH 

a  + 

-'f2U^X,  +  t,U^X2  +  U^ 

X3" 

j=l 
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L  ) 

* 
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A 
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- 
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V2 

Il-ui 
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A 

-'l'lX3*'|'3Xi  +  Ui'|'2X3-U^ 

X2^ 

j=l 

•J 

A 

J 

f 

A 

^2  _ 

'f3U,X2^  X$jUUv,*j 

m 

IS) 

+ 

'|',X2-'I'2Xi-U^ 

V2  _ 

X3*  X^jqj 

1=' 

j=i 

- 

- 

A 

J 

( 

A 

V2  _ 

V3U1X3+  Xf  jUl*v,*j 
j=l 

Ul 

J=l 

- 

- 

A 

J 

V2 

^  X  * 

+  U,  '|'2  X3  -  U,  '|'3  X2  +  2u,  Z  * 

j=1 
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-  U,  Xj  +  U,  '|'3  X,  )+ 


(  \ 

^  -.i 

Ui 

j=1 

- 

L  y 

* 

f  \ 

V2 

2ui  2<|>juVv,*j 

j=1 

^  + 

Ui 

V2 

j=l 

- 

- 

1  / 

^2  .2 

U,  'j',  X2  -  u,  'I'a  Xl  ^  2Ui  Z<t'jUi*v,*i))d3 

j=1 


5.  Panial  Velocities 

The  generalized  partial  velocities  and  angular  velocides  are  assembled  in  Table 
2.1  and  Table  2.2.  These  values  and  the  acceleration  equations  are  the  basis  of  the  formu¬ 
lation  for  the  generalized  inertia  forces. 


TABLE  2.1 

PARTIAL  VELOCITIES  OF  BODY  A  AND  BODY  B 
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TABLE  2.2 

PARTIAL  VELOCITIES  OF  BODY  C 


6.  Generalized  Inertia  Forces 

a.  Generalized  Inertia  Force  of  Body  A 

The  generalized  inertia  force  for  body  A  is  comprised  of  two  parts,  the 
first  is  due  to  the  moment  of  inertia  of  body  A’s  mass  center  about  the  83  axis.  The  second 
is  due  to  the  inertial  force  of  body  A’s  mass  center. 

A  /  naN  a**  na*< 

(Eqn.  2.25)  '  ("'^3  V  ^  •  2. 

where  J  3  is  the  moment  of  inertia  of  A*^  and  ITI  a  is  the  mass  of  body  A. 

A  A*^ 

Using  the  values  of  and  V  p  from  Table  2.1,  the  generalized  inertia 
force  of  body  A  is 

(Eqn.  2.26)  ^F,  =  J3  U ,  -  u  1 

There  is  only  one  equation  since  body  A  is  rigid  and  displays  only  rota¬ 
tional  motion. 

b .  Generalized  Inertia  Forces  of  Body  B 

Body  B  is  a  continuous  elastic  beam  with  a  mass  per  unit  length  of  p. 
Thus,  the  generalized  inertia  equations  are  an  integration  along  the  length  of  the  beam 

(Eqn.  2.27)  ®Fp  =  -  /^  (Yp^  pdx 

0 

The  result  is  an  equation  for  the  large  motion  and  one  equation  for  each 
deformable  degree  of  freedom  of  body  B. 

The  large  motion  equation  incorporates  only  first-order  terms,  while  the 


small  motion  equation  neglects  third-order  or  higher  terms. 


(Eqn.  2.28) 


=  -f  li,  '  P^P^2P, 


''^1  r  1  '^1 

v;.  L  xy  o  xx  .  o..  V 


^  Zuui  |Pi<t>rP2<t’ip  2u,  Zui*j  JPi<t>j+P2<|>il  pdx 


Pi*X<l>iqi  -<t>i  P2*X<t>?qi 


"  Xu'ui  |<l>j<t>^<^J<^^<t>j‘t>i  |-  ui  *  Pi<t>rP2<l>j"  X  qj 
i=l  i=l 


Vl 

>  2u,  £u,*j 

i=1 


Pdx 


c .  Generalized  Inertia  Forces  of  Body  C 

Body  C  was  modelled  with  lumped  masses  at  63  grid  points.  Therefore, 
the  generalized  inertia  equations  consist  of  a  summation  of  the  acceleration  dotted  with  the 
generalized  partial  velocity  times  the  mass 


(Eqn.  2.29) 


c  NX 
Yr-  ^ 


c 


There  are  I  +  v  j  +  V2  equations  of  motion  for  body  C.  Again,  the 
large  motion  equation  has  only  flrst-order  terms  and  the  subsequent  small  motion  equations 
neglect  third-order  and  higher  terms. 

(Eqn.  2.30) 
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(Eqn.  2.31) 
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'^2_ 
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lk(^iX2-'^'tx3 


■slktixsl  fui.vi-i  +k+rtkti 


The  generalized  active  forces  are  derived  from  the  forces  exerted  on  the  differ¬ 
ential  elements.  Only  the  internal  elastic  forces  and  the  external  forces  were  considered  in 
this  formulation.  Rotation  of  the  individual  beam  elements  was  neglected. 

The  equation  for  the  generalized  internal  forces  is  of  the  form 


(Eqn.  2.32) 


J  Yr®  ■  •  Gif'- 

0  X=0 


where  df  consists  of  internal  axial  and  shear  forces. 


The  shear  force  portion  of  df  f(X*  the  beam  is 


(Eqn.  2.33) 


shear 


dx  ^2 


where  y  = 

j=1 


The  axial  force  portion  of  df  for  the  beam  is 

a^u 

<^<^axlai  =  AE^dx  ai 


where  U  = 

j=l 


When  these  differential  forces  are  dotted  with  the  generalized  partial  velocity 
and  then  integrated  by  parts,  the  remaining  terms  are  simply  the  forces  due  to  strain  energy. 
The  generalized  internal  fences  are  reduced  to 
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(Eqn.  2.34) 


F,*j  =  -Bo).  rij  qj 


The  external  force  on  the  boom  was  the  torque  applied  about  83  axis  at  point  Q, 
the  fixed  end  of  the  boom.  This  torque  was  dotted  with  the  generalied  partial  angular 
velocity  to  produce 


(Eqn.  2.35) 


For  the  reflector  dish,  it  was  assumed  that  the  strain  energy  would  account  for 
the  internal  forces  and  no  external  fmces  were  applied  to  the  reflector. 


Thus,  the  generalized  acdve  forces  for  the  reflector  are 


(Eqn.  2.36) 


The  equations  of  Sections  B.6  and  B.7  were  gathered  to  obtain  the  following 


equations: 


(Eqn.  2.37) 
(Eqn.  2.38) 


(Eqn.  2.39) 


»  bf"  ♦  CfT- ,  =  0 

Bf".|  ♦  Bp".)  .  Cp".)  =  0  j  =  I . v, 

"“Fi-vi-k  ®F,*v,*k  =  0  k  =  1 . V2 


The  orthogonality  condition  for  the  boom  is: 


(Eqn.  2.40) 


A  i  I  ^  "  I"!  t  I  " 


1-2  r 

I  i: 

x=0  . 


and  for  the  reflector 


IIL  COMPUTER  IMPLEMENTATION 

A.  INTRODUCTION 

Computer  implementation  to  solve  these  equations  of  motion  involved  three  stages. 
Modal  analysis  was  conducted  using  the  NASA  Structural  Analysis  (NASTRAN). 
NASTRAN  is  a  general-purpose  digital  computer  program  for  the  analysis  of  various 
structures  [Ref.  7].  NASTRAN  provided  values  for  the  natural  frequencies,  mode  shapes, 
and  generalized  modal  masses  for  the  two  flexible  bodies. 

The  data  was  transferred  and  read  into  an  IBM/VS  FORTRAN  program.  Constants 
were  calculated  in  this  program,  which  then  were  inserted  into  the  initial  stage  of  a  Dynamic 
Simulation  Language  (DSL)  program.  DSL  is  a  digital  simulation  for  continuous  systems 
[Ref.  8].  DSL  was  used  to  solve  the  simultaneous  differential  equations. 

B .  MODAL  ANALYSIS 

A  multibody  dynamic  analysis  requires  special  consideration  of  the  boundary  condi¬ 
tions.  For  the  boom,  one  end  was  fixed  and  the  other  end  was  free  with  a  tip  mass 
attached.  The  tip  mass  was  equal  to  the  total  mass  of  the  reflector.  The  mode  shapes, 
natural  frequencies,  and  generalized  masses  were  calculated  using  17  grid  points  (Figure 
A.6).  Figures  A.7,  A.8,  A.9,  and  A.  10  depict  the  first  four  modes  of  the  boom. 

The  reflector  was  represented  by  63  grid  points  (Figure  A.  1).  Various  links  between 
grid  points  were  constrained  from  rotating  and  the  hinge  point  (Grid  63)  was  fixed. 
Figures  A.2,  A.3,  A.4,  and  A.5  depict  the  first  four  modes  of  the  reflection. 

Ten  modes  were  generated  by  NASTRAN  for  each  body.  A  preliminary  analysis 
between  two,  three,  four,  and  five  modes  for  each  body  was  conducted.  Four  modes  was 
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chosen  for  the  final  analysis  because  minimal  changes  occurred  between  the  four-  and  five¬ 
mode  cases. 


Modified  Given’s  method  (MGIV)  was  used  to  conduct  the  modal  analysis  of  the 
boom  and  the  reflector.  The  generalized  mass  was  normalized  to  equal  one  simplifying  the 
equations. 

The  natural  frequencies  for  the  boom  and  two  cases  of  the  reflector  are  tabulated  in 
Table  3.1.  In  the  results  section,  the  effects  of  reflector  flexibility  are  discussed.  The 
NASTRAN  programs  are  shown  in  Appendix  E. 

TABLE  3.1 

REAL  EIGENVALUES 


Mode 

Boom 

Natural 

Frequency 

(cycles) 

Reflector  1 
Natural 
Frequency 
(cycles) 

Reflector  2 

Natural 

Frequency 

(cycles) 

1 

.71674 

2.0937 

1.1635 

2 

.78059 

12.469 

6.9084 

3 

2.624 

15.929 

14.32 

4 

4.723 

15.929 

14.342 

C .  EVALUATION  OF  TIME  CONSTANTS 

Two  files  were  transferred  to  the  IBMA^S  computer.  Both  files  contained  ten  modes 
of  data  for  one  of  the  bodies.  All  ten  modes  are  read  into  a  FORTRAN  program,  but  the 
number  of  modes  used  to  calculate  constants  was  varied  by  a  parameter  statement  in  the 
beginning  of  the  program.  This  strategy  provided  easier  formatting  to  read  in  the  data. 


while  still  generalizing  the  program  and  making  the  computations  more  efficient.  The  pro¬ 
gram  is  shown  in  Appendix  E. 

The  equations  of  modon  were  systematically  evaluated  to  identify  terms  that  depended 
on  position  alone  and  were  repeated  several  times.  These  terms  were  isolated  from  the 
time-dependent  terms.  Subroutines  were  written  to  calculate  these  position  teims  and  were 
called  immediately  after  reading  in  the  mode  shape  data.  The  main  program  consisted  of 
multiplications  and  additions  between  system  constants  and  the  subroutine  outputs,  since 
the  subroutines  carried  out  the  required  position  integrations  or  summations. 

The  reason  for  having  separate  programs  is  that  all  the  position  calculations  can  be 
completed  beforehand.  Thus,  the  DSL  program  need  only  be  concerned  with  time  step 
calculations.  The  number  of  constants  and  their  dimensions  depended  on  the  DSL  imple¬ 
mentation  and  will  be  discussed  in  that  section. 

D.  DSL  IMPLEMENTATION 

With  the  separate  FORTRAN  program  to  evaluate  the  time  constants,  the  DSL  pro¬ 
gram  is  relatively  simple.  The  equations  of  motion  are  put  in  the  form 

[A]  ii  =  a 

where  [A]  is  a  matrix  of  the  coefficients  of  the  acceleration  terms  (the  left-hand  side  of  the 
equations  of  motion),  while  ^  is  a  vector  of  all  other  terms  (the  right-hand  side  of  the 
equations  of  motion). 

These  coefficients  were  supplied  by  the  FORTRAN  program.  The  A  matrix  and  B 
vector  were  assembled  in  the  derivative  section  of  the  DSL  program.  Two  subroutines 
from  LINPACK  were  called  to  decompose  the  A  matrix  and  solve  the  matrix  equation. 
DGEFA  uses  gaussian  elimination  to  decompose  A  and  DGESL  solves  the  matrix 
equation. 
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DSL  provided  many  alternatives  to  solve  a  time-stepping  problem  for  simultaneous 
differential  equations.  A  Runge-Kutta  fifth-order  method  with  variable  step  size  was 
selected  because  it  was  self-starting,  stable,  and  accurate.  Although  other  methods  may 
have  been  more  efficient  computationally,  accuracy  was  a  primary  selection  criterion.  The 
program  is  shown  in  Appendix  E. 

The  results  were  printed  out  in  a  file  and  plotted  using  the  TEK618.  Plotting  was 
accomplished  via  the  GRAFAEL  command,  which  is  a  particularly  powerful  accessory  of 
the  DSL  program. 


IV.  RESULTS 


Three  parameters  of  the  LFMR  system  were  varied  in  the  computer  simulation  to 
investigate  the  effects  of  these  parameters  on  pointing  error  of  the  LFMR  system  during  a 
spin-up  procedure.  Deflections  and  slope  changes  at  the  boom  tip  (grid  point  63)  and  grid 
points  1  and  4  of  the  reflector  were  compared  for  this  purpose.  The  locations  of  these  grid 
points  are  shown  in  Figure  4.1.  The  first  analysis  compared  the  effects  of  changing  the 
reflector  flexibility,  while  the  boom  properties  remained  the  same.  The  second  analysis 
compared  vibrational  amplitudes  of  three  cases  with  different  in-plane  orientation  of  the 
reflector.  These  orientations  were  accomplished  by  lining  up  the  a  i  ,02,33  coordinate 
system  with  the  C 1 ,02,03  coordinate  system  and  rotating  the  C 1  and  O3  axes  about  the  82 
axis  to  the  desired  angle.  The  third  analysis  compared  the  case  of  a  particular  in-plane 
orientation  (-155*  case)  with  the  case  of  the  same  in-plane  orientation  plus  slight  out-of- 
plane  tilt  (+5*).  The  out-of-plane  orientation  was  achieved  by  rotating  the  reflector  about 
the  03  axis. 

The  same  spin-up  procedure  was  used  in  all  computer  simulation  runs  and  the  torque 
applied  to  the  LFMR  system  during  the  spin-up  procedure  is  shown  in  Figure  4.2. 

A .  MATERIAL  PROPERTIES  AND  PARAMETERS 

The  material  used  to  model  the  boom  and  reflector  is  Isotropic  Graphite-Epoxy 
Composite  (T3()0/5208  (0/90/457-45)3)  [Ref.  9].  The  inertia  property  of  the  electronics  box 

was  obtained,  assuming  it  to  be  a  uniform  two-foot  cubic  body.  The  properties  of  the 
boom  and  reflectors  are  presented  in  Table  4.1.  Reflector  1  is  the  baseline  design  of  the 
reflector  used  throughout  the  analyses.  Reflector  2  is  a  more  flexible  model  of  the  reflector 
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and  was  used  in  the  reflector  flexibility  analysis.  The  geometric  parameters  for  the  boom 
and  two  models  of  the  reflectors  are  given  in  Table  4.2. 


TABLE  4.1 

BOOM  AND  REFLECTOR  PROPERTIES 


Property 

Boom 

Reflector  1 

Reflector  2 

Modulus  of  elasticity 
(E) 

10.1E  +  6psi 

10.  IE  +  6  psi 

10.  IE  +  6  psi 

Poisson’s  Ratio 

(V) 

.25 

.25 

.25 

x-sectional  area 
(A) 

1.1868  in” 

.1425  in2 

.1425  in 

Inner  diameter 

Rod  element 
(di) 

2.5  in 

1.04859  in 

.8177  in 

Outer  diameter 
rod  element 
(di) 

3.0  in 

1  in 

.75  in 

Thickness  of  rod  element 

.5  in 

.04859  in 

.0677  in 

1.522  E-4  slugs/in^ 

Area  moment 

of  inertia 
(A) 

1.2231  in^ 

.0165  in^ 

.00853  in^ 

TABLE  4.2 


GEOMETRIC  PARAMETERS 


Parameter 

1  ] :  length  of  lower  boom 

162  in 

I2:  length  of  upper  boom 

126  in 

cx:  angle  between  1  \  and  83  axis 

34.5* 

angle  between  1  ^  and  I2 

98.76* 

M^:  weight  of  electronics  box 

501b 

M^:  weight  of  reflector 

37.5  lb 

B  .  EFFECTS  OF  THE  REFLECTOR  FLEXIBILITY 

Two  cases  of  reflector  flexibility  were  chosen  for  comparison.  As  shown  in  Table 
4.1,  Reflector  1  was  stiffer  than  Reflector  2.  The  difference  in  the  angular  velocity  is 
shown  in  Figures  B.l  and  B.2.  Greater  amplitude  oscillations  occur  after  the  applied 
torque  is  removed  for  the  more  flexible  reflector.  This  effect  is  due  to  larger  vibrations  of 
the  more  flexible  reflector  (Reflector  2). 

The  fourth  generalized  coordinate  (Q<^(t))  of  the  boom  was  chosen  to  represent  the 
reflector  flexibility  effects  on  the  boom’s  generalized  coordinates.  Figures  B.3  and  B.4 
depict  this  boom  generalized  coordinate  for  the  two  cases.  They  show  a  definite  dynamic 
interaction  between  the  boom  and  reflector.  The  more  flexible  reflector  displayed  greater 
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fluctuations  in  the  boom  generalized  coordinate,  which  translated  into  greater  vibrations  at 
the  boom  tip.  The  vertical  deflections  of  the  boom  tip  are  shown  in  Figures  B.7  and  B.8. 

The  third  reflector  generalized  coordinate  (03(1))  was  chosen  to  represent  the  effects 
of  reflector  flexibility  on  the  reflector’s  generalized  coordinates.  Comparison  of  Figures 
B.5  and  B.6  shows  a  negative  shift  in  the  graph  of  the  more  flexible  reflector.  The  ampli¬ 
tude  of  the  reflector  generalized  coordinate  is  unchanged. 

The  rotational  vibrations  about  the  vertical  axis  are  compared  in  Figures  B.9  and 

B.  10.  The  more  flexible  reflector  caused  greater  amplitude  rotational  vibration  at  the  boom 
tip. 

C.  EFFECTS  OF  IN-PLANE  ORIENTATION  CHANGES 

As  mentioned  previously,  cases  of  different  in-plane  orientations  between  the  reflec¬ 
tor  and  the  boom  were  analyzed:  Three  cases  of  angle  rotations  about  the  82  axis  were 
compared:  -135*,  -145*,  and  -155*.  The  angular  displacement  for  these  three  orientations 
is  shown  in  Figures  C.1-3.  The  angular  displacements  oscillated  after  the  applied  torque 
was  set  to  zero.  An  increasing  amplitude  trend  in  angular  displacements  was  apparent  as 
the  orientation  was  a  less  negative  angle. 

The  horizontal  deflection  at  the  boom  tip  represents  the  effects  of  reflector  orientation 
on  the  boom.  Figures  C.4-6  show  that  greater  amplitude  vibration  occurred  for  the  case  of 
-135*  reflector  orientation  in  the  horizontal  direction. 

The  rotations  at  the  boom  tip  about  the  horizontal  axis  are  shown  in  Figures  C.7-9. 
Greater  fluctuations  occurred  for  the  -135*  orientation.  The  larger  amplitude  results  for  the 
case  of  -135*  orientation  is  due  to  the  location  of  the  center  of  mass  of  the  reflector.  The 
radius  of  rotation  is  larger  for  the  -135’  orientation.  Thus,  the  moment  of  inertia  is 
increased. 
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The  reflector  orientation  also  affects  the  reflector  vibration.  Figures  3.10-12  show 
the  horizontal  displacement  of  grid  point  1.  Higher  frequency  components  appear  in  the 
-155*  case,  but  the  greater  amplitude  vibrations  occur  in  the  -135*  case. 

D.  EFFECT  OF  OUT-OF-PLANE  ORIENTATION  OF  REFLECTOR 

The  case  of  -155*  in-plane  cmentaQtxi  of  the  reflector  is  compared  to  the  case  with  the 
same  155*  in-plane  rotation  plus  a  5*  tilt  out  of  plane.  The  first  boom  generalized  coordi¬ 
nate  (Q  ^(t))  was  the  only  boom  coordinate  affected  by  the  tilt  out  of  plane.  The  tilt  pro¬ 
duced  a  positive  shift  in  (Q  j  ( t))  value  as  shown  in  Figures  D.  1  and  D.2.  The  oscillations 
were  same  in  amplitude. 

The  fourth  reflector  generalized  coordinate  (Q^Ct))  was  chosen  to  represent  the 
effect  of  the  out-of-plane  orientation  on  the  reflector  generalized  coordinates.  Figures  D.3 
and  D.4  show  that  vibrations  with  larger  amplitude  appeared  for  the  out-of-plane 
orientation. 

The  out-of-plane  orientation  deflections  for  the  boom  and  reflector  were  the  most 
affected  by  the  out-of-plane  tilt  of  the  reflector.  Figures  D.5  and  D.6  show  the  effects  at 
the  boom  tip.  The  flrst  boom  generalized  coordinate  was  the  predominant  component 
resulting  in  a  positive  shift  in  the  out-of-plane  displacement  of  the  boom.  Figures  D.7  and 
D.8  show  the  effects  at  grid  point  4  with  the  first  boom  generalized  coordinate  producing  a 
positive  shift.  The  reflector  generalized  coordinates  produced  a  larger  amplitude  vibration 
for  the  out-of-plane  rotation. 
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A.  CONCLUSIONS 


The  purpose  of  this  research  was  to  develop  a  realistic  model  of  the  LFMR  sensor 
system.  This  goal  was  successfully  met  by  applying  Kane’s  method  to  a  multibody  model. 
Efficient  computer  simulation  was  achieved  The  programs  were  general  enough  to  use  up 
to  ten  modes  for  each  body,  but  the  results  indicate  that  superimposition  of  four  modes  for 
each  body  was  sufficient. 

The  flexibility  of  the  reflector  did  affect  the  vibration  of  the  boom  and  was  related  to 
pointing  error.  The  multibody  system  allowed  the  effects  of  reflector  flexibility  to  be 
investigated  without  reanalyzing  the  entire  system. 

The  program  also  applied  easily  to  different  orientations  of  the  reflector.  Modal  anal¬ 
ysis  did  not  have  to  be  repeated  as  in  the  case  of  a  single-body  model. 

B .  RECOMMENDATIONS 

Application  of  this  research  could  include  a  study  of  the  effects  of  a  full  range  of  dif¬ 
ferent  configurations  on  the  system  dynamics.  Also,  damping  effects  can  easily  be  incor¬ 
porated  in  this  work  to  study  the  dynamic  interaction  of  each  subsystem.  Extension  of  this 
research  could  include  a  study  of  the  effects  of  a  flexible  body  connected  to  the  spacecraft’s 
main  body  on  three-dimensional  motion  and  control. 


APPENDIX  A 


MODAL  ANALYSIS  OF  REFLECTOR  AND  BOOM 


Figures  A.l  through  A.  10  show  the  undeformed  as  well  as  the  first  four  mode  shapes  of  the 


reflector  and  boom.  The  figures  are  on  pages  84-93. 
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APPENDIX  B 


REPRESENTATIVE  EFFECT  OF  REFLECTOR  FLEXIBILITY 

Figures  B.l  through  B.IO  show  representative  effects  of  two  different  reflector  flexibilities 
on  the  dynamic  response  of  the  system.  The  figures  are  on  pages  96-105. 
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APPENDIX  C 


REPRESENTATIVE  EFFECT  OF  IN-PLANE  REFLECTOR  ROTATION 

Figures  C.l  through  C.12  show  representative  effects  of  three  different  in-planc  orientations 
of  the  reflector  on  the  dynamic  response  of  the  system.  The  figures  are  on  pages  106-117. 
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APPENDIX  D 

REPRESENTATIVE  EFFECT  OF  OUT-OF-PLANE  REFLECTOR  ROTATION 


Figures  D.l  through  D.8  show  representative  effects  of  an  out-of-plane  orientation  of  the 
reflector  cmi  the  dynamic  response  of  the  system.  The  figures  are  on  pages  1 18-125. 
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APPENDIX  E 


PROGRAMS 


id  nastran  Natalie 
sol  3 
time  500 
cend 

title=Nomal  Mode  Analysis  of  LFMR  Boom 

subtitle=Research  on  N-ROSS  DYNAMICS 

method® 10 

displacement=all 

spc=100 

output(plot) 

set  1  =  all 

origin  5 , 0 . ,0. 

axes  y,2,x 

view  -30. ,20. ,0. 

ptitle=Undeformed  Boom 

find  scale, origin  5, set  1 

plot  set  1, label  grid 

ptitle=Mode  Shape 

maximum  deformation  15 

find  scale, origin  5, set  1 

plot  modal, deformation, 0,1  THRU  10, set  l,pen  4, shape 
begin  bulk 
$ 

SEIGR, Sid, METHOD, fl,f 2, ne,nd, , ,+ 

$  +,N0RM,g,c 
$ 

eigr,10,MGIV.0.0,100.0, ,10 
$ 

$CORD2C,cid,rid,al ,a2,a3,bl,b2,b3,+ 

S  +,cl,c2,c3 
$ 

cord2r,l, ,0. ,0. ,0, ,0. ,0. ,1. ,+2r 
+2r , 1 . ,0, , 1 . 

$ 

$GRID, id,cp,xl ,x2,x3,cd,ps,seid 
$ 

grid, 1 , 1 , 0 . ,0 . ,0 . 

=,*(!),*, *(10. 1953), *.*(14. 83427),==  $ 

*7 

grid, 10, 1,91. 75781 ,0. ,133.5084 
=,*(!),», *(-13. 1085), =,*(12. 3356),==  $ 

=6 

grid, 18, 1,0. ,0. ,100. 

$ 

$CBEAM , eid , pid , ga , gb , gO/xl , x2 , x3 
$ 

cbeam,l,3,l,2,18 
=,*(!),=, *(!),*(!),==  $ 

^14 
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$ 

$PBEAM,pid,midl ,A, II , 12 , J ,nsm 
$ 

pbeam, 3 ,1,1 .1868,1.2231 , 1 .2231 , ,2.4462 

$I1AT1  ,mid,E,G,Nu,Rho,A,Tref  ,Ge 
$ 

$  A  =  Thermal  Expan.  Coef. 

$  Tref  =  Ref.  Temp. 

$  Ge  =  Damping  Coef.  (=  2  x  C/Co  ) 

$ 

matl,l,101.+5, ,0.25,1522.-7 
$ 

$CONM2,eid,g,cid,m,xl,x2,x3, ,+ 

$  +,111,121,131,132,133 
$ 

conm2, 50, 17. 1.9. 4502-2 
$ 

spcl, 100, 123456, 1,18 
$ 

enddata 


id  nastran  Natalie 
sol  3 
time  500 
cend 

title=Nomal  Mode  Analysis  of  LFMR  Reflector 

subtitle=Research  on  N-ROSS  DYNAMICS 

method=10 

displacement*all 

spc=l00 

output(plot) 

set  1  *  all 

origin  5,0. ,0. 

axes  z,x,y 

view  34.27,-23.17,0. 

ptitle=Undeformed  Reflector 

find  scale, origin  5, set  1 

plot  set  1, label  grid 

ptitle=Mode  Shape 

maximum  deformation  15 

find  scale, origin  5, set  1 

plot  modal. deformation,!  THRU  10, set  l,pen  4, shape 
begin  bulk 
$ 

$EIGR, Sid, METHOD, fl,f2,ne,nd, . ,+ 

S  ♦,NORM,g,c 
$ 

eigr .10, MGIV, 0.0,100.0, ,10 
$ 

$CORD2C,cid.rid,al,82,a3,bl,b2,b3,+ 

S  •t-,cl,c2,c3 
S 

ccrd2r,:,,0.,0.,0,,0.,0.,l.,+2r 
+2r.l. ,0. ,1. 

cord2c ,2.1,0.,0.,0.,0.,0.,1. ,+ab 
fab, 10. ,0, ,10. 

$ 

SGRID, id,cp,xl ,x2,x3,cd,ps,seid 
$ 

grid, 1,2, 116. 81, 15, ,40.83 
*,’*(1),  =  ,  =  ,*{30.),»=  $ 

=  10 

grid, 13, 2, 89. 26, 15. ,36.71 
=,*(!),», *.*(30.),==  $ 

=  10 

grid, 25, 2, 60. 86, 15. ,32.45 
=,*(1) ,=,=,*(30. ) ,==  $ 

=  10 

grid, 37, 2, 31. 5, 15. ,30.13 
=  , *(!),  =  ,  =  .*(30.),==  $ 

=10 

grid, 49, 2, 63. 8, 15. ,4.72 
=,*(1),«.=,*(30.),=»  $ 

=10 

.grid, 61, 1,0.  ,0.  ,31.67 


grid, 62, 1,0. ,0. ,27.65 
grid, 63, 1,0. ,0. ,0. 
grid, 64, 1 ,0 . ,0 . , 50. 
grid, 65, 2, 50. ,0. ,0. 

$ 

$CBEAM , eid , pid , ga , gb , gO/xl , x2 , x3 
$ 

cbeam, 1 ,1,1,13,64 

=,"(1) ,=,*(1) .*(1) .==  $ 

=10 

cbeam, 13, 1,13, 25, 64 

=,*(1),=,*(1), *(!),==  $ 

=10 

cbeam ,25,1,25,37,64 

=,*(1),=,*(1), *(!),==  $ 

=10 

cbeam, 37, 1,37, 62, 64 
=,*(!),=, *(!),==  $ 

=10 

cbeam, 49, 2, 63, 62, 65 
cbeam, 50, 2, 62, 6 1,6 5 
$ 

crod,l,10,25,49 

=  ,*(!),  =  , Ml), *(!),==  $ 

=10 

crod,13,ll,l,49 

=,*(1),=,*(1), *(!),==  $ 

=10 

crod,25,ll,25,63 
=,*(!),=, *(!),==  $ 

=10 

crod,37,ll,49,63 
=,*(!),=, *(!),==  $ 

=  10 

crod,49,ll,l,2 
=  ,*(!),  =  , Ml), *(!),==  $ 

=9 

crod, 60 ,11,12,1 
crod,61,ll,13,14 

=  ,*(!),  =  ,*(!), Ml), ==  $ 

=9 

crod, 72, 11, 24, 13 
crod, 73, 11, 25, 26 
=  ,*(!),  =  ,*(!), Ml), ==  $ 

=9 

crod, 84, 11, 36, 25 
crod, 85, 11, 37, 38 

=  ,M1),=,M1),M1),==  $ 

=9 

crod, 96, 11, 48, 37 
crod, 97, 11, 25, 50 
=  ,*(!),  =  ,*(!), Ml), ==  $ 

=9 

.crod, 108, 11, 36, 49 
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crod,109,ll,49,26 
=,*(1).=,*(1), *(!),==  $ 

=9 

crod, 120, 11 ,60,25 
$ 

$PBEAI1,pid,midl  ,A,I1 ,12,J,nsm 
$ 

pbeam,] , . .0,1452 ,0.01647 .0 .01647 , ,0.03295 
pbeap-..2,l  3 . 927 , 3 . 1907 ,3 . 1907  .  ,6 .3314 
prcd.lvj,!  ,0.03  236,1 .882-3,1.0 
prod, 11 ,1,0.01094,1.911-5,1.0 
$ 

SMATl ,mid.E,G,Nu,Rho,A,Tref ,Ge 
$ 

S  A  =  Thermal  Expan.  Coef. 

$  Tref  =  Ref,  Temp. 

$  Ge  =  Damping  Coef.  (=  2  x  C/Co  ) 

$ 

matl, 1,101. +5, ,0.25,1522.-7 
$ 

$C0NM2,eid,g,cid,m,xl ,x2,x3, ,+ 

S  +,111,121,131,132,133 
$ 

conm2, 1,1, 2, 0.4562-3 

=,*(1). *(!),==  $ 

=10 

conm2, 13, 13, 2, 0.5432-3 
=  , Ml), *{!).==  $ 

=10 

co.Tm2, 25, 25, 2, 0.4122-3 

=,*(1), *(!),==  $ 

=10 

conm2, 37 ,37,2,0.4438-3 
=  , Ml), *(!).==  $ 

=10 

conm2, 49, 62,1 ,0,626-2 
conm2 ,50 ,63 ,1 ,0 ,649-2 
$ 

spcl , 100 , 123456 ,63  thru  65 
$ 

enddata 
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THIS  PROGRAM  IS  DESIGN  TO  READ  IM  MODAL  INFORMATION  FROM 
A  NASTRAN  DATA  FILE  AND  CALCULATE  TIME  CONSTANTS  FOR  A  DSL 
PROGRAM.  TWO  FILES  ARE  READ  CONTAINING  THE  DATA  FOR  THE  LFMR 
BOOM  AND  REFLECTOR.  A  PARAMETER  STATEMENT  IS  USED  TO  DEFINE 
THE  NUMBER  OF  GRID  POINTS  AND  MODES  FOR  BOTH  BODIES  IN  THE 
BEGINNING  FOR  A  MORE  GENERAL  PROGRAM. 

THE  FOLLOWING  PARAMETERS  ARE  USED  IN  THIS  PROGRAM  : 

NPT.NPTC  :  THE  NUMBER  OF  GRID  POINTS  FOR  BODY  B  AND  C 
RESPECTIVELY 

M,MC  ;  THE  NUMBER  OF  MODES  FOR  BODY  B  AND  C 

RESPECTIVELY  TO  BE  USED  TO  CALCULATE  CONSTANTS 


J3 

B2MA 

P1,P2,P3 

X1,X2,X3 

AIIASS 

CMASS 

BMASS 

RHO 

DX 

FIBX,FIBY,FIBZ 

FICX.FICY.FICZ 

PSI1,PSI2,PSI3 

OHGAB , OMGAC 

P1P2IN 

PSFSIN 

BINII 

PHBXY 

PHZS 

CMS 

XISMS 

X2SMS 

X3SMS 

XICM 

X2CM 

X3CM 

X12MS 

X23HS 

X13MS 

PBX,PBY,P5Z 

PBXS,PBYS,PBZS 

PBXY 


:  THE  MOMENT  OF  INERTIA  OF  BODY  A  ABOUT  THE 
SPIN  AXIS 

:  THE  INERTIA  DUE  TO  THE  DISPLACEMENT  OF 
BODY  A'S  CENTER  OF  MASS  FROM  THE  SPIN  AXIS 
:  THE  COORDINATES  OF  THE  GRID  POINTS  ON 
BODY  B 

!  THE  COORDINATES  OF  THE  GRID  POINTS  ON 
BODY  C 

:  THE  TOTAL  MASS  OF  BODY  A 
:  AN  ARRAY  OF  THE  MASS  AT  EACH  GRID  POINT  OF 
BODY  C 

A  DUMMY  ARRAY  TO  READ  IN  MODAL  DATA 
THE  HASS  PER  UNIT  LENGTH  OF  BODY  B 
THE  INTERVAL  DISTANCE  BETWEEN  GRID  POINTS  ON 
BODY  B 

:  MATRICES  OF  BODY  B'S 
SHAPES  (ROWS  ARE  DOFS, 

:  Mi^TRICES  OF  BODY  B'S 
SHAPES  (SAME  FORMAT  AS 


INTEGRAL 

INTEGRAL 

INTEGRAL 

INTEGRAL 


OF 

OF 

OF 

OF 


SUMMATION  OF 
SUMMATION  OF 
SUMMATION  OF 
SUMMATION  OF 
SUMMATION  OF 
SUMMATION  OF 
SUMMATION  OF 
SUIUIATIOM  OF 
SUMMATION  OF 
SUMMATION  OF 
SUMMATION  OF 


* 

★ 

* 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 


TRANSLATIONAL  MODE 
COLUMNS  ARE  MODES) 
TRANSLATIONAL  MODE 
FIBX.FIBY.FIBZ  ) 

:  ARRAYS  OF  THE  ROTATIONAL  MODES  OF  THE  HINGE 
POINT  BETWEEN  BODY  B  AND  BODY  C 

THE  AlIGULAR  FREQUENCY  OF  BODIES  B  AND  C 
INTEGRAL  OF  Pl**2+P2**2’'-RHO*DX 

(P1*FIBX  +  P2*FIBy)*RHO*DX 
(Pl^FIBY  +  P2*FIBX)*RHO*DX 
(FIBXd,  J)*FIBY(I,K)*RHO*DX 
FIBZ(I, J)*FIBZ(I ,K)*RHO*DX 
GRID  POINT  HASSES  OF  BODY  C 
Xl**2  *  CMASS 
X2**2  *  CMASS 
X3**2  *  CMASS 
XI  *  CMASS 

X2  *  CMASS  * 

X3  *  CMASS  * 

X1*X2*CMASS  * 

X2*X3*CMASS  * 

X1*X3*CHASS  * 

FICX,Y,Z  *  CMASS  * 

SUmiATON  OF  FICX,Y,Z(I,J)*FICX,Y,Z(I,K)*CMASS* 
SUMIIATION  OF  FICX(I ,  J)*FICY(I  ,K)*CMASS  * 


* 

A 

* 

* 

* 

* 

* 

5^ 

* 

* 

* 

A: 

* 

* 

* 

:k 

•k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 


X1FBX,Y,Z 

X2PBX.Y,Z 

X3PBX,Y,Z 

UlDll 

U1D12 

U1D13 

UJDll 

U1D21 

U1D22 

UKDll 

U1D31 

U1D32 

UJD21 

UKD21 

B12UUJ 

B12UUK 

BUZ  21 
BU222 

BU223 


B22UUJ 

B22UUK 

B3U211 

B3U213 

B32UUJ 

B32UUK 


SUMMATION  OF  X1*FICX, Y,2*CMASS  * 
SUMMATION  OF  X2*FICX,Y,Z*CMASS  * 
SUMMATION  OF  X3*FICX, Y,Z*CMASS  * 
A  (1,1)  CONSTANT  TERM  * 
A  (1.1)  COEFFICIENTS  OF  BODY  B  GENERALIZED  * 
COORDINATES  * 
A (1,1)  COEFFICIENTS  OF  BODY  C  GENERALIZED  * 
COORDINATES  * 
A(1,J+1)  * 
A(l+J,l)  COEFFICIENTS  OF  BODY  B  GENERALIZED  * 
COORDINATES  * 
A{l+J,l)  COEFFICIENTS  OF  BODY  C  GENERALIZED  * 
COORDINATES  * 
A(l,l+M+K)  * 


A(l+M+K,l)COEFFICIENTS  OF  BODY  B  GENERALIZED  * 
COORDINATES  * 

A(l+M+K.l)COEFFICIENTS  OF  BODY  C  GENERALIZED 
A(l+I,l+J)  * 

A(l+J,l+M+K)  « 

B(l)  COEFFICIENTS  OF  2U1  TIMES  BODY  B'S  * 

GENERALIZED  SPEEDS  * 

B(l)  COEFFICIENTS  OF  2*U1  TIMES  BODY  C'S  * 
GENERALIZED  SPEEDS  * 

B(l+J)  Ul**2  CONSTANT  TERM  * 

B(l+J)  COEFFICIENTS  OF  Ul**2  TIMES  BODY  B'S  * 
GENERALIZED  COORDINATES  * 

B(l+J)  COEFFICIENTS  OF  Ul**2  TIMES  BODY  C'S  * 
GENERALIZED  COORDINATES  * 

ALSO  B(l+M+K)  COEFFICIENTS  OF  Ul’''*2  TIMES  * 
BODY  B'S  GENERALIZED  COORDINATES  * 

B(l+J)  COEFFICIENTS  OF  2*U1  TIMES  BODY  B'S  * 
GENERALIZED  SPEEDS  * 

B(l+J)  COEFFICIENTS  OF  2*U1  TIMES  BODY  C'S  * 
GENERALIZED  SPEEDS  * 

B(l+M+K)  CONSTANT  COEFFICIENTS  OF  Ul**2  * 

B(l+M+K)COEFFICIENTS  OF  Ul**2  TIMES  BODY  C'S  * 
GENERALIZED  COORDINATES  * 

B(l+M+K)COEFFICIENTS  OF  2*U1  TIMES  BODY  B'S  * 
GENERALIZED  SPEEDS  * 

B(l+M+K)COEFFICIENTS  OF  2*U1  TIMES  BODY  C'S  * 
GENERALIZED  SPEEDS  * 


IMPLICIT  REAL*8(A-H,0-Z) 

INTEGER  NPT,NPTC,M,MC,I,J,K 
PARAMETER  (NPT=17 ,NPTC=63  ,M=10,MC=10) 

REAL*8  J3 

DIMENSION  Pl(NFT) ,P2(NPT) ,P3(NPT) ,PISQ(NPT) ,P2SQ(NPT) ,RPFMIX(NPT ,M 
+) ,P1SP2S(NPT) ,RPSSQ(NPT) ,CHASS(NPTC) ,X1(NPTC) , 

+X2(MPTC) ,X3(NPTC) ,FIBZ(NPT,H) ,PHI2C(NPTC) , 

+FI3X(NPT,M) ,FIBY(NPT,H) ,FICX(NPTC,MC) ,FICY(MPTC,MC) ,FICZ(NPTC 
+,MC) ,P1FBX(NPT,M) ,P2N(NPT) ,P2FXN(NPT,M) ,P1FY(NPT,M) ,PFIMIX(NPT,M> 
+,P2FBY(NPT,M) ,P1P2FA(NPT,M) ,RPSFA(NPT,M) ,PSFSIN(H) ,PSI1 (M) , PSI2 (H) 
+,PSI3(M) ,FBX(MC) ,PBy(MC) ,X1PBX(MC) ,X2PBY(HC) ,X2PBX(MC) ,X1PBY(HC) , 
+X1PB2(MC) ,X2PBZ(MC) ,X3PBX(MC) ,X3PBY(MC) ,X3PBZ(MC) , 


+U1D12(M) ,U1D13(MC) ,UJD11(M) .UKDll(MC) ,B12UUJ(M) ,B12UUK(HC) , 
+BU221(M),B3U211(HC).U1D21(M,M),U1D22(M,MC),UJD21(M,H),UKD21(M,MC), 
+BU222(M,M) ,BU223(M,MC) ,B22UUJ(M,M) ,B22UUK(M,MC) ,B32UUJ(M,HC) , 
+U1D31(H,MC) ,B3U213(MC.MC) .B32UUK(MC,MC) , 

+PHIPHI(NPT) ,PHBXY(M,M),PHZS(M,M) ,PBZ(MC) ,PBXS(MC,MC) , 

+PBYS(HC,HC) ,PBXy(HC,MC),BINlI(ll) ,U1D32(MC,HC) ,PBZS(MC,MC) , 
+FREQB(H) ,OMGAB(M) ,FREQC(MC) .OMGAC(HC) ,BMASS(MPT) 

DATA  J3, B, AMASS, DX,RHO/12.0, 12. , .1295,18. ,1.806310-4/ 

READ  IN  CONSTANTS 

CALL  MOD.JiTB  (M . MPT .  PI , P2 .  P3 ,  FIBX ,  FIBY ,  FIBZ ,  PSI 1 ,  PSI2 ,  PSI3 ,  BMASS , 
+FREOB,OMGAB) 

CALL  MODATC (MC , NPTC , XI , X2 , X3 , FICX , FICY , FICZ , CMASS , 

+FREQC,OMGAC) 

CALL  ALL  CONSTAiMT  SUBROUTINES 

CALL  BINT11(NPT.RH0.DX.P1,P2.P1SQ,P2SQ,P1SP2S,RPSSQ,P1P2IN) 

CALL  BI.MT12  (NPT, M, PI  ,P2, FIBX, FIBY, P1FBX,P2FBY,P1P2FA,RPSFA, 
+RHO,DX,PSFSIM) 

CALL  BINTII  (NPT, M, PI ,P2, FIBX, FIBY, RH0,DX,P2N,P2FXN,P1FY, 
+FFIMIX,RPF!IIX,BIN1I) 

CALL  CALPXY (NPT , M , F IBX , FIBY , RHO , DX , PHIPHI , PHBXY) 

CALL  FIBZSQ (NPT , M , FIBZ . PHIPHI , DX , RHO , FHZS ) 

CALL  SUIIASS (NPTC, CMASS, CMS) 

CALL  XSQC!1S  (NPTC  ,  CMASS  ,  XI  .XI ,  XI SMS ) 

CALL  XSQCMS  (MPIC, CMASS,. X2,X2,X2SMS) 

CALL  XSQCMS  ( NPTC, CMASS ,X3,X3 ,X3SHS) 

CALL  SUMC  (NPTC, XI, CMASS, XICM) 

CALL  SUMC  (NPTC,X2,C.MASS,X2CM) 

CALL  SUMC  (NPTC, X3, CMASS, X3CM) 

C.ALL  XSQCMS  (NPTC .CMASS , XI  ,X2 , XI 2MS ) 

CALL  XSQCMS  (NPTC, CMASS, X2,X3,X23MS) 

CALL  XSQCMS  (NPTC , CMASS , XI , X3 , XI 3MS ) 

C.ALL  C.ALFBM  (NPTC, MC, FICX, CMASS, PBX) 

CALL  CALFBM  (NPTC, MC, FICY, C^tASS,PBY) 

CALL  CALPBM  (NPTC, MC, FICZ. CMASS, FBZ) 

CALL  P5SYMS  (NPTC, MC, FICX, CMASS, FHI2C,PBXS) 

CALL  FBSYMS  (NPTC, MC, FICY, CMASS, PHIZC.FBYS) 

CALL  PBSYMS  (NPTC, MC, FICZ, CMASS, PHI2C,PBZS) 

CALL  CALFSS  (NPTC ,HC , FICX , FICY , CMASS , PHI2C ,PBXY) 

CALL  CALXPB  (NPTC ,MC , XI , FICX, CMASS , XIPBX) 

CALL  CALXPB  (NPTC ,MC , XI , FICY , CMASS , XIPBY) 

CALL  CALXPB  (NPTC ,HC , XI , FICZ, CMASS ,X1PBZ) 

CALL  CALXPB  (NPTC ,MC , X2 , FICX, CMASS , X2PBX) 

CALL  CALXPB  (NPTC, MC ,X2 , FICY, CMASS , XZPBY) 

CALL  CALXPB  (NPTC ,MC,X2 , FICZ , CMASS ,X2FBZ) 

CALL  CALXPB  (NPTC, MC,X3, FICX, CMASS, X3PBX) 

CALL  CALXPB  (NPTC, MC,X3, FICY, CMASS, X3PBY) 

CALL  CALXPB  (NPTC, MC ,X3 , FICZ, CMASS ,X3PBZ) 

B2MA  =  B’^B’^AMASS 

H1H2CH  =  (P1(NPT)*P1(NPT)  +  P2(NPT)^P2(NPT) )*  CMS 
UlDl 1  =J3+B2MA+PlP2IN+HlH2CM+XlSMS+X2SMS+2 . * (PI (NPT ) *X1CM+P2 (NPT ) 
+  *X2CM) 


FXTER  =  FIBX(NPT,J)*(P1(NPT)*CMS  +  XICM) 

FYTER  =  FIBY(NPT,J)*(P2(NPT)*CMS  +  X2CM) 

SIITER  =  -PSI1(J)*(P2(HFT)*X3CK  +  X23HS> 

SI2TER  =  PSI2(J)*(P1(NFT)*X3CM  +  X13MS) 

SI3TER  =  PSI3(J)*(P2(MPT>*X1CM  -  PI (NPT)*X2CM) 

BU221(J)  =  PSFSIM(J)+FXTER+FYTER+SI2TER+SI3TER+SI1TER 
U1D12(J)  =  2.  *  BU221(J) 

10  CONTINUE 

DO  11  K  =  1,MC 

U1D13(K)=2.*(P1(NPT)*PBX(K)  +X1PBX(K)  +  P2 (NPT)*PBY(K)  +  X2PBY{K)) 

11  CONTINUE 

DO  12  J  =  l.M 

FXTER  =  FIBX(NPT,J)*(P2(NPT)*CMS  +  X2CM) 

FYTER  =  FIBY(NPT.J)*(P1(NPT)*CMS  +  XICM) 

SIITER  =  PSI1(J)*(P1(NPT)*X3CM  +  X13MS) 

SI2TER  =  PSI2(J)*(P2{NPT)*X3CH  +  X23MS) 

SI3TER  =  PSI3(J)*(P1(NPT)*X1CM+X1SMS+X2SMS+P2(NPT)*X2CM) 
UJDll(J)  =  BIN1I(J)-FXTER+FYTER-SI1TER-SI2TER+SI3TER 

12  CONTINUE 

*CALCULATES  A (1,3)  AND  PART  OF  A (3.1) 

DO  13  K  =  1,MC 

UKDll(K)  =  -PBX(K)*P2(NPT)-X2PBX(K)+PBY(K)*P1(NPT)  +  XIPBY(K) 

13  CONTINUE 

DO  14  J  =  l.M 

FXTER  =  FIBX(NPT.J)*(P1(NPT)*CMS  +  XICM) 

FYTER  =  FIBY(NPT,J)*(P2(NFT)*CMS  +  X2CM) 

SIITER  =  PSI1(J)*(P2(NPT)*X3CM  +  X23MS) 

SI2TER  =  PSI2(J)*(P1(NPT)*X3CM  +  X13MS) 

SI3TER  =  PSI3(J)*(P1(NPT)*X2CM  +  P2(NPT)*X1CM) 
B12UUJ(J)=-PSFSIN(J)-FXTER-FYTER-SIITER+SI2TER-SI3TER 

14  CONTINUE 

DO  15  K  =  1,MC 

B12UUK(K)=-P1(NPT)*PBX(K)-X1PBX(K)-P2(NPT)*PBY(K)-X2PBY(K) 
B3U211(K)  =  -B12UUK(K) 

15  CONTINUE 
*  1+J  EQUATION 

A 

DO  16  I  =  l.M 
DO  17  J=1,M 

F.XTER=FIBX(NPT,J)*(PSI1(I)*X3CM-PSI3(I)*X1CM-FIBY(NPT,I)*CHS) 
FYTER=FIBY(NPT,J)*(PSI2(I)*X3CM-PSI3(I)*X2CM+FIBX(NPT,I)*CIIS) 
SI1TER=PSI1(J)*(  (PSI3(I)*P2(NPT)  +  FIBX(MPT , I ) )*X3CM  + 
+(PSI1(I)*P1(NFT)  +  PSI2(I)*P2(NPT))*X2CM  +  PSI2(I)*(X3SMS+X2SHS) 
++PSIl(I)*Xi2!:S) 

SI2TER=PSI2(J)*(  (PSI3(I)*P1(NPT)  -  FIBY(NPT, I ) ) *X3CM  + 
+(PSI1(I)*P1(NPT)  +  PSI2(I)*P2(NPT))*X1CM  +  PSIl (I ) *(X3SMS+X1SMS ) 
++PSI2(I)*X12MS) 

SI3TER=PSI3(J)*(  (-PSI3(I)*P1(NPT)  +  FIBY(NPT , I ) )*X2CM  + 
+(PSI3(I)’"P2(MPT)+FIBX(NPT,I))*X1CH-PSI1(I)*X23MS+PSI2(I)*X13MS) 
U1D21(J,I)=PHBXY(I, J)-PHBXY(J.I)+FXTER+FYTER-SI1TER+SI2TER+SI3TER 
17  CONTINUE 

16  CONTINUE 
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DO  18  J  =1  ,M 
DO  19  K  =  1,MC 

FXTER=FIBX(NPT,J)*PBY(K) 

FYTER=FIBY(NPT,J)*PBX(K) 

SI1TER=-PSI1(J)*(X3P3X(K)+P1(NPT)*PBZ(K)+X1PBZ(K)) 

S I 2TER=-PS I 2 ( J ) * ( X3PBY( K ) +P2 ( HPT ) *PBZ ( K ) +X2PBZ ( K ) ) 
SI3TER=PSI3(J)*(2.*(X2PBY(K)+XiPBX(K))+P2(NPT)*PBY(K)+Pl(NPT)* 
+  PBX(K)) 

U1D22(J,K)=-FXTER  +  FYTER  +  SIITER  +  SI2TER  +  SI3TER 

19  CONTINUE 
18  CONTINUE 

DO  20  1=  1,M 
DO  21  J  =  I,M 

FXTER=FIBX(NPT,J)*(PSI2(I)*X3CM  -  PSI3(I)*X2CH) 

+  +FIBX(NPT,I)*(PSI2{J)*X3CM  -  PSI3(J)*X2CM) 

FYTER=FIBY(NPT, J)*(-PSI1(I)*X3CM  +  PSI3 ( I ) *X1CM) 

+  +FIBY(NPT,I)*(-PSI1(J)*X3CM  +  PSI3{J)*X1CM) 

FZTER=FIBZ(NPT, J)*(PSI1(I)*X2CM  -  PSI2 ( I ) *X1CM) 

+  +FIBZ(NPT,I)*(PSI1(J)*X2CM  -  PSI2{J)*X1CM) 

SI1TER=PSI1(J)*(PSI1(1)*(X3SMS+X2SHS)-PSI2(I)*X12MS) 

+  +PSI1(I)*{PSI1(J)*(X3SMS+X2SMS)-PSI2{J)*X12MS) 

SI2TER=FSI2(J)*(PSI2(I)*(X3SMS+X1SMS)-PSI3(I)*X23HS) 

+  +PSI2(I)*(PSI2(J)*(X3SMS+X1SMS)-PSI3(J)*X23MS) 

SI3TER=PSI3(J)*(PSI3(I)*(X2SMS+X1SMS)-PSI1(I)*X13HS) 

+  +PSI3(I)*(PSI3(J)*(X2SMS+X1SMS)-PSI1(J)*X13MS) 

IF(I  .EQ.  J)THEN 

UJD21(I ,J)=1 .+FXTER+FYTER+FZTER+SI1TER+SI2TER+SI3TER 
ELSE 

UJD21(I,J)=FXTER+FYTER+F2TER+SI1TER+SI2TER+SI3TER 

UJD21(J,I)=UJD21(I,J) 

END  IF 

21  CONTINUE 

20  CONTINUE 

*  CALCULATES  A(2,3)  AND  A(3,2) 

DO  22  K  =  1,MC 
DO  23  J  =  1  ,11 

FXTER=FIBX(NPT, J)*PBX(K) 

FYTER=FIBY(NPT,J)*PBY(K) 

FZTER=FIBZ(NPT, J)*PBZ(K) 

SI1TER=FSI1(J)*(X2PBZ(K)-X3PBY(K)) 
SI2TER=PSI2(J)*(X3PBX(K)-X1PBZ(K)) 
SI3TER=PSI3(J)*(X1FBY(K)-X2?BX(K) ) 

UKD21 ( J , K)=  FXTER+FYTER+FZTER+SI1TER+SI2TER+SI3TER 
23  CONTINUE 

22  CONTINUE 

* 

DO  24  J  =  1,M 
DO  25  I  =  1,M 

FXTER=FIBX(NPT, J)*(-PSI3(I)*X2CH+PSI2(I)*X3CM) 
FYTER=FIBY(NPT,J)*(-PSI1(I)*X3CM+P3I3(I)*X1CH) 
SI1TER=PSI1{J)*(-FIBY(NPT,I)*X3CM+PSI1(I)*(X3SMS-X2SMS'P2(NPT) 
+*X2CI1)+PSI2(I)*(P1(NPT)*X2CM+X12HS)+FSI3(I)*P1(NPT)*X3CH) 

SI2TER=PSI2(J)*(FIBX(MPT,I)*X3CN+PSI2(I)*(X3SHS-X1SHS-P1(NPT) 


59 


+*XlCH)+PSIl(I)*(P2(NFT)*XlCH+Xi2MS)+PSI3(I)*P2(NPT)*X3CH) 

SI3TER-PS13(  J)*(-FIBX(NFT,  I)  ‘•X2CIl-FSI3(I)*(r2(MPT)*>:2CMfPl(NPT) 
+*X1CM)-PSI1(I)*(X13MS)-PSI2(I)*X23MS+FIBY(NPT,I)*X1CM) 

IF  (I  .EQ.  J)  THEN 

BU222 { J , I ) =1 . -PHZS ( J , I ) -CMS*FIBZ (NPT , I ) *FIBZ ( NPT , J ) +FXTER+FYTER 
+  +SIiTER+SI2TER+SI3TER 

ELSE 

BU222 ( J , I )=-PHZS { J , I ) -CHS*FIBZ(MPT , I ) *FIBZ (NPT , J ) +FXTER+FYTER 
+  +SI1TER+SI2TER+SI3TER 

END  IF 

25  CONTINUE 

24  CONTINUE 


DO  26  K=1,MC 
DO  27  J  =  1,M 
FXTER=FIBX(NPT , J) *PBX(K) 

FYTER=FIBY(NPT,J)*FBY(K) 

SI1TER=-PSI1(J)*(X3PBY(K)+P2(NPT)*PBZ(K)+X2PBZ(K)) 
SI2TER=PSI2(J)*(X3PBX(K)+P1(NPT)*PBZ(K)+X1PBZ(K)) 
SI3TER=PSI3(J)*(P2(NPT)*PBX'K)-P1(NPT)*PBY(K)) 
BU223(J,K)=FXTER+FYTER+SHTER+SI2TER+SI3TER 
27  CONTINUE 

26  CONTINUE 


DO  28  J=1,H 
DO  29  1=1, M 

FXTER=FIBX(NPT,J)*(FIBY(NPT,I)*CMS+PSI3(I)*X1CM-PSI1(I)*X3CM) 
FYTER=-FIBY(NPT,J)*(FIBX(NPT,I)*C11S+PSI2(I)*X3CM-PSI3(I)*X2CM) 
FZTER=FIBZ(NPT,J)*(.5*PSI3(I)*X3CM) 
SIlTER=PSIl(J)*{FIBX(NPT,I)*X3CM+PSI2(I)*X3SMS-PSI3(I)*1.5*X23nS) 
SI2TER=PSI2(J)*(FIBY(NPT,I)*X3Cn-PSIl(I)*X3SMS+PSI3(I)*1.5*X13HS) 
SI3TER=PSI3(J)*(-FIBY(NPT,I)*X2CM-PSI2(I)*X13MS+PSI1(I)*X23MS- 
+  FIBX(NPT,I)*X1CM) 

B22UUJ(J,I)=PHBXY(I,J)-FHBXY(J,I)  +  FXTER  +  FYTER  +  FZTER  + 

+  SI1TER+SI2TER+SI3TER 

29  CONTINUE 
28  CONTINUE 

DO  30  K  =  1,MC 
DO  31  J  =  1,M 
FXTER=  FIBX(NPT, J)*PBY(K) 

FYTER=  -FIBY(NPT, J)*PBX(K) 

SI1TER=  PSI1(J)*X3PBX(K) 

SI1TRK=  SIITER  -FSI 1 ( J ) * (PI (NPT) *PBZ(K)+X1PBZ(K) ) 

SI2TER=  PSI2(J)*X3PBY(K) 

SI2TRK=  SI2TER  -PSI2( J)*(P2(NPT)*PBZ(K)+X2PBZ(K) ) 

SI3TER=  -PSI3(J)*(X1PBX(K)+X2PBY(K) ) 

SI3TRK  =  PSI3(J)*(P1(NPT)*PBX(K)+P2(NPT)*PBY(K)) 

B22UUK(J,K)  =  FXTER  +  FYTER  +  SIITER  +  SI2TER  +  SI3TER 
B32UUJ(J,K)  =-B22UUK(J,K)+.5*X3PBZ(K)*PSI3(J) 

U1D31(J,K)  =  FXTER  +  FYTER  +  SIITRK  +  SI2TRK  +  SI3TRK 
31  CONTINUE 

30  CONTINUE 

*  l+N+K  FOU;*TIGM  (URDU  IS  IRST  TERM;  U1D31  IS  2HD;IDENTITY  IS  3RD) 


GO 


*  B3  {B3U2H  IS  IRST  TERH;BU223  IS  U2  Q(I)  TERM; 

*  B32UUJ  IS  2U  *U(1+I)TERM 

DO  32  K  =  1,MC 
DO  33  J  =K,HC 

IF  (K  .EQ.J)  THEN 

B3U213  (J,K)  =  1.  -  PBZS{J,K) 

ELSE 

B3U213  (J,K)  =  -PBZS(J,K) 

B3U213(K,J)  =  B3U213(J,K) 

END  IF 

33  CONTINUE 

32  CONTINUE 

DO  34  K  =  1,MC 
DO  35  J  =  1,MC 

B32UUK(J,K)  =  PBXY(K,J)-PBXY{J,K) 

U1D32(J,K)  =  -  B32UUK(J,K) 

35  CONTINUE 

34  CONTINUE 

C  PRINT  OUT  RESULTS  IN  FILE  25 

WRITE  (25, ' (lOX, ' 'AllC ' ) ' ) 

WRITE  (25,100)  UlDll 
100  FORMAT  (D15.6) 

WRITE(25, ' (lOX, ' 'AllQ' ' ,9X, ' 'AJIC ' ,8X, ' ‘BIQ' ’ ,9X, ' 'BJUIC  ' ,9X, 

+ ' ' OMGAB ' ' ) ' ) 

DO  1  J=  1,M 

17RITE(25,200)J,U1D12(J) ,UJD11(J) ,B12UUJ(J) ,BU221(J) ,OMGAB(J) 

1  CONTINUE 

200  FORMAT  (I4,5D13.6) 

NRITE(25, ' (lOX, ' 'BIQB' ' ,9X, ' 'AllQB' ' ,8X, ' 'AIKC ' ,9X, ' 'BKUIC ' ,9X, 

+ ' ' OMGAC ' ' ) ' ) 

DO  2  K  =  1,MC 

HRITE(25,200)K,B12UUK(K) ,U1D13(K) ,UKD11 (K) ,B3U211(K) ,OMGAC(K) 

2  CONTINUE 

WP.ITE(25,  '  (14X,  '  'AJJC*  MIX ,  ” AJQ "  ,  12X ,  '  'BJUIQ'  ’  ,13X,  '  'BJ2UQ’  ' )  '  ) 
DO  3  I  =  1,M 
DO  4  J  =  i,n 

WRITE(25,500)J,I,UJD21(J,I) ,U1D21(J,I) ,BU222(J, I) ,B22UUJ(J,I) 

4  CONTINUE 

3  CONTINUE 

500  FORMAT  (2I4,4D15.6) 

WRITE (25, ' (12X, ' 'AKIQB' ' ,9X, ' 'BKUIQB' ' ,9X, ' ’BK2UQB' ' ) ' ) 

DO  5  K  =  1,MC 
DO  6  I  =  1,MC 

WP.ITE(25,300)I,K,U1D32(I,K)  ,B3U213(I,K)  ,B32UUK(I,K) 

6  CONTINUE 

5  CONTINUE 

300  FORnAT(2I4,3D15.6) 

WRITE(25, ' (12X, ' 'AJIQB' ' ,9X, ' 'AKJC' ' ,12X, ' 'BKUIQ' ' ,11X, ' 'BJ2UQB' ' ) 
+  ') 

DO  7  K  =1,MC 
DO  3  J  =  1,H 

WRITE (25, 500 ) J, K, U1D22 ( J,K) ,UKD21(J,K) ,BU223(J,K) ,B22UUK(J,K) 

8  CONTINUE 

7  CONTINUE 
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WRITE(25,  '  (12X,  '  'AKIQ'  '  ,9M,  '  •BI<2UQ'  '  )  '  ) 

DO  98  K  =  1,MC 
DO  99  J  =  1,M 

WRITE  (25,400)  J.K,U1D31(J,K) ,B32UUJ(J,K) 

99  CONTINUE 

93  CONTINUE 

400  FORMAT  (2I4,2D15.6) 

DO  94  J  =  1,H 
DO  95  I  =1,NPT 

WRITE(25,900)FIBX(I,J) ,FIBY(I,J) ,FIBZ(I,J) 

95  CONTINUE 

94  CONTINUE 

DO  96  J  =  1,MC 
DO  97  I  =  l.NPTC 

WRITE (25, 900)  FICX( I , J) , FICY( I , J) , FICZ( I , J) 

97  CONTINUE 

96  CONTINUE 

DO  93  J  =  1,M 

WRITE  (25,900)  PSIl ( J) ,PSI2 ( J) , PSI3 ( J) 

93  CONTINUE 

DO  92  1=  l.NPTC 

WRITE  (25,900)  XI ( I ) , X2 ( I ) , X3 ( I ) 

92  CONTINUE 

900  FORMAT  (3D20.6) 

END 

•k  k 

^  VECADD  * 

*  (  ADD  TWO  ARRAYS  TO  PRODUCE  AN  ARRAY  )  * 

k  k 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

SUBROUTINE  VECADD  (MPT, A, B , ARC) 

REAL--8  A(NPT)  ,B(NPT)  ,ARG(NPT) 

DO  10  I  =  1,NPT 

ARG(I)  =  A(I)  +  B(I) 

10  CONTINUE 

RETURN 
END 

*  MULTSV  ^ 

*  (  MULTIPLIES  A  SCALAR  TINES  A  VECTOR  )  * 

SUBROUTINE  MULTSV  (NPT , A ,B , ARG) 

RE AL*8  A , B ( NPT ) , ARG ( NPT ) 

DO  10  I  =  1,NPT 

ARG(I)  =  A  *  B(I) 

10  CONTINUE 

RETURN 
END 

*  MATADD  * 

*  (  ADDS  TWO  MATRICES  TO  PRODUCE  A  MATRIX  )  * 


Sl'BP.OUTIHE  IIATADD  ( MPT , H ,  A ,  B  .  AP.G ) 

REALMS  A(!IPT,M)  ,B(MPT,tl)  ,ARG(IIPT,M) 

DO  10  J  =  1,M 

DO  20  I  =  1,NPT 

AP.Od.J)  =  A(I,J)  + 

20  CONTIMUE 

10  CONTIMUE 

RETURN 
END 

*  SUBROUTINE  I1ULT2S  * 

^  (  CALCULATES  A  COLN  TIMES  A  ROW  VECTOR  )  * 

xAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA^cAAAA 

SUBROUTINE  MULT2S  (  MPT,A,B,ARG) 

REALMS  A(NPT) ,B(MPT) ,ARG(NPT) 

BO  10  I  =  1,NPT 

ARG{I)  =  A(I)  *  E(I) 

10  CONTINUE 

RETURN 
END 

AAAAAAAkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

*  MULTSM  * 

*  (  MULTIPLIES  A  SCALAR  TIMES  A  MATRIX  )  * 

AAAAAA,tAkAAAAAAAAAAAAAAAAAAAAAAAAA:^AAAAAAAAAAAAAAAAAAAAAkAAAAAAAAAAAAAAA 

SUBROUTINE  MULTSM  (NPT,M,A,B,ARG) 

REALMS  A,B{NPT,M) ,ARG(MPT,M) 

DO  10  J  =  1,M 

DO  20  I  =  1,NPT 

ARG(I,J)  =  A*  B(I,J) 

20  CONTINUE 

10  CONTINUE 

RETURN 
END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAs^^^^t*'**^*' 

*  SUBROUTINE  MULTVM  * 

^  (  MULTIPLIES  VECTOR  AND  MATRIX  TO  PRODUCE  A  MATRIX  )  * 

AA-rAAAA  VAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA  +  AAAAAAAAAAAAAAAAA 

SUBROUTINE  NULTVH  (NPT ,H, A , B , ARG ) 

REAL'^S  AUIPT)  ,B(NPI,M>  ,ARG(NPT,n) 

DO  10  I  =  1,NPT 
DO  20  J  =1,N 
ARG(I.J)  =  A(I)  *  B(I,J) 

20  CONTIMUE 

10  CONTINUE 

RETURN 
END 

AAAAAArA-i-AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

*  SUBROUTINE  FIE2SQ  * 

(  CALCULATES  THE  INTEGRAL  OF  PIB2  ( I  ,  J )  *PIBZ  ( I ,  K )  *RHO*DX  )  * 

7:AAAAA--AA':>AAAA:AAAAAAAAAAAAAAAAAA»AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE  FIBZSQ (NPT . N  .  FIEZ ,  FHIFHI ,  D-” .  RHO  ,  PHZS  ) 

FEAL*8  FIBZ(NFT.H)  .FHIFHKNPT)  ,PHZS(M,N)  ,RHO,DX 
DO  10  K  =  1  ,11 
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DO  20  J  =  K,M 
DO  30  I  -  l.NFT 

PHIPHI(I)  =  FIBZ{I ,J)*FIBZ(I ,K)*RH0 
30  CONTIMUE 

CALL  TRAF2{NPT,DX,PHIPHI,SUH) 

PHZS(J,K)  =  SUM 
PHZS(K, J)=PHZS(J,K) 

20  CONTINUE 

10  CONTIMUE 
RETURN 
END 

*  TRAPZ  * 

*  (  PERFORMS  INTEGRATION  APPROXIMATION  USING  TRAPEZOIDAL  RULE)  * 

SUBROUTINE  TRAPZ (NPT , DX . F , INT ) 

EEAL*3  DX,F(NFT) ,INT,SUM 
SUM  =  0. 

DO  10  1=  1,8 

IF(I  .EQ.  1  .OR.  I  .EQ.  NPT)  THEN 
FACTOR  =  1. 

ELSE 

FACTOR  =  2. 

EIIDIF 

SUM  =  SUM  +  FACTOR*F(I) 

10  CONTINUE 

HIT  =  DX/2.  *  SUM 

RETURN 

END 

^  MATTRP  * 


*  (INTEGRATION  APPROXIMATION  OF  A  MATRIX  USING  TRAPEZOIDAL  RULE)  * 

SUBROUTINE  MATTRP (NPT , M , DX , F , INT ) 

REAL-^S  DX,F(NPT,M)  ,INT(H)  .SUM 
DO  5  J  =  1,2 
SUM  =  0. 

DO  10  1=  1,8 

IF(I  .EQ.  1  .OR.  I  .EQ.  S)  THEN 
FACTOR  =  1. 

ELSE 

FACTOR  =  2. 

END  IF 

SUM  =  SUM  +  FACTOR*F(I,J) 

CONTINUE 

INT(J)  =  DX/2.  *  SUM 
CONTINUE 
RETURN 
END 

*  SUMC  * 

*  (  PERFORMS  SUMMATION  OF  GRID  POINT  HASSES  OF  BODY  C  TIMES  * 

*  AN  ARBITRARY  VECTOR  OF  EQUAL  DIMENSION  )  * 


SUBROUTINE  SUI1C(NPTC ,  F,  CHASS  ,  SUM) 

REALMS  F(MPTC) ,aiASS(MPTC) , SUN 
SUM  =  0. 

DO  10  I  =  1,NPTC 

SUM  =  SUM  +  F(I)*CMASS(I) 

10  CONTINUE 
RETURN 
END 

*  HAT SUM  * 

'  (  PERFORMS  SUMMATION  OF  THE  GRID  POINT  MASES  OF  BODY  C  TIMES  * 

*  AN  ARBITRARY  MATRIX  OF  EQUAL  DOF)  * 

^■k^A-k'^-k-ki'-k'k'k-k-k-k'^ri^-k'kiiiKi^iK-ki^'k-f'.-kTfi-k-ki^A-k-kK-k'k-k-k'kitiK-kiKiKii'k-k-k'k-kTV'kiK'k'kiK'k^iv'kiK^'k'ki^'k-k-k^ 

SUBROUTINE  MATSUM  (MPT , M , F , CHASS , SUM ) 

REAL*8  F(NPT,H) ,CMASS(MFT) ,SUH(M) 

DO  10  J=  1,M 

DO  20  I  =  1,NPT 

SUH(J)  =  SUM(J)  +  F(I, J)*CMASS(I) 

20  CONTINUE 

10  CONTINUE 

RETURN 
END 

BINTll  * 

*  (  PERFORMS  INTEGRATION  OF  (Pl**2  +  P2**2*RHO*DX)  * 

SUBROUTINE  BINTll (HPT ,RHO , DX , PI , P2 , PISQ , P2SQ , P1SF2S , RPSSQ , P1P2IH) 
REAL*3  Pl(NPT) ,P2(NPT) ,PiP2IN,PlSQ(NPT) ,P2SQ(NPT) ,RPSSQ(MPT) 

REALMS  P1SP2S(NPT) ,FHO,DX 
CALL  HULT2S(NPT,P1 ,P1,P1SQ) 

CALL  HULT2S(NPT,P2,P2,P2SQ) 

CALL  VECADD (NPT , PISQ , P2SQ , P1SP2S ) 

CALL  MULTSV(NPT,RH0,P1SP2S,RFSSQ) 

CALL  TRAPZ(NPT,DX,RPSSQ,P1P2IN) 

RETURN 

END 

*  BINT12  * 

*  (  PERFORMS  INTEGRATION  OF  (P1*FIBX  +  P2*FIBY) *RHO*DX  )  * 

SUBROUTINE  BINT12  (NPT ,M , PI , P2 , FIBX , FIBY , PIFBX , P2FBY , P1F2FA , RPSFA , 
+RHO,DX,PSFSIN) 

REAL*8  Pl(NPT) , P2 (NPT) , FIBX(NPT , M) , FIBY(NPT , M) , PIFBX (NPT , H) , 
+P2FBY(NPT,M) ,P1P2FA(NPT,H) ,PSFSIN(M) ,RPSFA(NPT,H) ,RHO,DX 
CALL  MULTVM(NPT, PI , FIBX, PIFBX) 

CALL  MULTVM(NPT , P2 , FIBY , P2FBY ) 

CALL  MATADD (NPT, PIFBX, P2FBY,P1P2FA) 

CALL  i:ULTSM(NFT,H,RHO,PlP2FA, RPSFA) 

CALL  MATTRP ( NPT , M , DX , RPSFA , PS FS IN ) 

RETURN 

END 

BINTll 

*  (  PERFORMS  INTEGRATION  OF  (Pl-^FIBY  -  P2*FIBX)  *RHO*DX  ) 
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'-/rf  UV  . 


fit;*, 


SUBROUTINE  BINTil  (NPT  ,H ,  PI  ,  P2  ,  FI''X  ,  FIB'/ ,  RHO  ,  Da  ,  P2N  ,  P2FXH  ,  FI  FT’ , 
+FrIMIX,RPFMIX,BINII) 

REALMS  PI(NPT) ,P2(NPT) ,FIBX(HFT,H) ,FIBY(NPT,M) ,P2H(NPT) ,BINII(M) 

+  ,?2FXN(NPT,M)  ,P1FY(MPT,I1)  ,PFII1IX(NPT,H)  ,RPFHIX(NPT,1-I)  ,RHO,DX 
DO  10  I  =  I, NPT 
P2M(I)  =  -P2(I) 

10  CONTINUE 

CALL  MULTVM  (NPT ,H , P2N , FIBX , P2FXN) 

CALL  NULTVH  (NPT ,M , PI , FIBY , PIFY) 

CALL  MATADD  (NPT ,M , P2FXN , PIFY , PFIMIX) 

CALL  IIULTSH  (NPT ,M , RHO , PFINIX , RPFHIX) 

CALL  HATTRP  (NPT,M,DX,RPFMIX,BIN1I ) 

RETURN 

END 

^  CALPXY 

INTEGRAL  OF  (  FIBX(  I ,  J) *FIBY(  I ,  K)  *RHO*DX)  * 

AAAAAAAAAAAAAAAAAAAnAAAnAnAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE  CALPXY  (  NPT , M , FIBX , FIBY , RHO , DX , PHIPHI , PHBXY) 

REAL*8  FIBX(NPT,M) ,FIBY(NPT,M) ,PHBXY(M,M) ,PHIPHI(NPT) , RHO, SUM 
+  ,DX 

DO  10  K  =  1,M 
DO  20  J  =  1,M 
DO  30  I  =  1,NPT 

PHIFHI(I)  =  FIBX(I,J)’"FIBY(I,K)*RHO 
30  CONTINUE 

CALL  TRAPZ (NPT, DX, PHIPHI, SUM) 

PHBXY (J,K)  =  SUM 
20  CONTINUE 

10  CONTINUE 
RETURN 
END 

*  SUMASS  * 

*  (SUMMATION  OF  THE  GRID  POINT  MASSES  OF  BODY  C  )  * 

SUBROUTINE  SUMASS  (MPTC , CMASS , CMS ) 

REALMS  CMASS (NPTC) , CHS 
CMS  =  0. 

DO  10  I  =1,NPTC 

CMS  =  CHS  +  CMASS (I) 

10  CONTINUE 

RETURN 
END 

AAAAAAAAAAAAAAAAAAAAxAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAnA 

*  XSQCMS 

*  (  SUMMATION  OF  THE  MULTIPLICATION  OF  TWO  ARBITRARY  X'S  TIMES  * 

*  THE  GRID  POINT  MASSES  OF  BODY  C  )  ■" 

SUBROUTINE  XSQCMS  (NPTC , CMASS , XA , XB , XSHS ) 

F.EAL-3  XA(NPTC)  ,XB(NPTC)  ,  CMASS (iNPTC )  ,XSMS 
ZSMS  =  0. 

DO  10  I  =I,NPTC 


XSHS  =  XSHS  +  XA(I)*XB(I)*C!1A5S(I) 

10  CONTINUE 
RETURN 
END 

*  CALPBH  * 

*  (SUMMATION  OF  AN  ARBITRARY  MODE  SHAPE  MATRIX  TIMES  THE  GRID  * 

*  POINT  MASSES  OF  BODY  C  )  * 

SUBROUTINE  CALPBH  (NPTC,MC,FIC,CHASS ,PB) 

REALMS  FIC(NPTC,MC)  ,CI1ASS(NPTC)  ,PB(MC)  ,SUM 
DO  10  J  =  i,nc 
SUM  =  0. 

DO  20  I  =  1,NPTC 

SUM  =  SUM  +  FIC(I,J)*CHASS(I) 

20  CONTINUE 

PB(J)  =  SUM 
10  CONTINUE 

RETURN 
END 

^  CALPBS  * 

*  {  SUMMATION  OF  THE  MULTIPLICATION  OF  TWO  ARBITRARY  MODE  SHAPES  * 

TIMES  THE  GRID  POINT  HASSES  OF  BODY  C  )  * 

SUBROUTINE  CALPBS  (NFTC .MC , FICA , FICB , CHASS , PHI2C , PBS ) 

REALMS  FICA(NPTC,MC) ,FICB(MPTC,MC) ,CMASS(NPTC) ,PBS(MC,MC) 

REAL*3  PHI2C(NPTC) ,SUH 
DO  10  K  =  1,MC 
DO  20  J  =  1,MC 
DO  30  1  =  1,NPTC 
PHI2C(I)  =  FICA(I,J)*FICB(I,K) 

30  CONTINUE 

CALL  SUHC (NFTC , PHI2C , CHASS , SUM ) 

F3S  (J,K)  =  SUM 
20  CONTINUE 

10  CONTINUE 
RETURN 
END 

FBSYMS  * 

^  (SUMMATION  A  MODE  SHAPE  MATRIX  TRANSPOSED  TIMES  THE  MODE  SHAPE  * 
MATRIX  TIMES  THE  GRID  POINT  MASSES  OF  BODY  C)-  A  SYMMETRIC  MATRIX  * 

SUBROUTINE  PBSYMS  (NPTC ,MC , FICA , CMASS , PHI2C , PBS ) 

REAL"^8  FICA(NFTC,HC)  , CHASS (NPTC)  , PBS  (MC , MC) 

REAL-3  PHI2C(NPTC) ,SUn 
BO  10  K  =  1,MC 
DO  20  J  =  K,MC 
DO  30  I  =  i,NPTC 
PHI2C(I)  =  FICA(I,J)^FICA(I,K) 

30  CONTINUE 

CALL  SUMC ( NPTC, PHI2C, CHASS , SUM) 

FBS  (J,K)  =  SUM 


FBS  (K,J)  =  PBS  (J.K) 

20  CONTINUE 

10  CONTINUE 
RETURN 
END 

*  CALMFB  * 

*  (  SUNNATION  OF  ANY  X*ANY  MODE  SHAPE  MATRIX*THE  GRID  POINT  HASSES  * 

*  OF  BODY  C  )  * 

KAAAAKAjirx:TA>:AAAAAAAAAAAAAAAAAAAAAAAAAAAKAAAAAAA^AAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE  CALXPB  (NPTC , MC , X , FIC , CHASS , XFB ) 

REAL’^8  FIC (NPTC, MC) ,CMASS (NPTC) ,XPB(HC) ,SUM,X(NPTC) 

DO  10  J  =  1,KC 
SUM  =  0. 

DO  20  I  =  l.NPTC 

SUM  =  X(I)*FIC(I,J)*CMASS(I)  +  SUM 
20  CONTINUE 

XPB(J)  =  SUM 
10  CONTINUE 

RETURN 
END 

AAAAAAA^'VkkAAAAAAkAkAAAAAAAAAAAAAAAAAAAAAAAAAAAxAAAAAAAA 

*  IlODATB:  READ  AND  TAYLOR  MODAL  DATA  OF  BODY  B  * 

AAAAAAAAAAAAAAAAAkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE  HODATB  (M ,NPT , XI , X2 , X3 , FIBX , FIBY , FIBZ , 

+  PSl , PS2 , PS3 , CMASS , FREQ , OMEGA) 

C 

IMPLICIT  REAL*3(A-H,0-Z) 

DIMENSION  FREQ(M) ,ONEGA(M) ,CMASS(NPT) 

DIMEN'SIOM  XKI.'PT)  .X2(NPT)  ,X3(NPT) 

DIMENSION  FIBX(NFT,M) ,FIBY(NFT,M) ,FIBZ(NPT,M) 

DIMENSION  FSl(M) ,PS2(M) ,PS3(M) 

DItlEMSION  FSI1(100,30) ,PSI2(100,30) ,PSI3(100,30) 

DIMENSION  VB(3) ,VC(3) ,VO(3) , VN(3) , C(3 , 3 ) 

C 

REWIND  8 
REWIND  21 
C 

PI  =  4,0  *  ATAN(l.DO) 

TPI  =  2.0  *  PI 
FRAD  =  PI/180.0 
C 

READ (8, 100)  ICOR,IMAS,ICMOD,ICTR,NGRID,NMOD 
C 

IF(ICTR  .LE.  1)  GO  TO  3 
IF(ICTR  .GE,  3)  GO  TO  2 
READ(8,104)  ALP 
GO  TO  3 

2  READ(8,106)  VB(1) ,VB(2) ,VB(3) , 

+  VC(1) ,VC(2) ,VC(3) 

C 

3  IF(IMAS  .EQ.  0)  GOTO  9 
C 


DO  5  I=1,II!AS 


n  o  o  n  o 


READ(3,150)  J1,J2,BMAS 
DO  4  J=J1,J2 

4  CMASS(J)  =  BIIAS 

5  CONTINUE 
C 

9  DO  10  I=1,NGRID 

10  READ(a,200)  X1(I),X2(I),X3(I) 

C 

IF(ICOR  .EQ.  0)  GO  TO  30 
C 

DO  20  I=1,NGRID 
R  =  X1{I) 

THTA  =  X2(I)  *  FRAD 
X1(I)  =  R  *  COS{THTA) 

20  X2(I)  =  R  *  SIN(THTA) 

C 

30  DO  50  I=l,NHOD 

READ (8, 300)  IMOD 
READ (8, 400)  FREQ (IMOD) 

OMEGA (IMOD)  =  TPI  *  FREQ (IMOD) 

C 

DO  40  J=1,NGRID 

READ(8,500)  FIBX(J,IMOD) ,FIBY(J, IMOD) ,FIBZ(J, IMOD) , 
+  PSIl(J,IMOD),PSI2(J,IMOD) ,PSI3(J,IMOD) 

40  CONTINUE 

FSi(I)  =  FSI1(NPT,I) 

FS2(I)  =  FSI2(N?T,I) 

PS3(I)  =  PSI3(N?T,I) 

50  CONTINUE 

IF(ICTR  ,LE.  1)  GO  TO  80 
IFdCTR  .GE.  3)  GO  TO  75 

2- DIMENSIONAL  COORDNATE  TRANSFORM 

ALP  =  ALP  *  FRAD 
CP  =  COS (ALP) 

SP  =  SIN(ALP) 

DO  60  I=1,NGRID 
X  =  XI (I) 

Z  =  X3(I) 

XI (I)  =  X  *  CP  +  Z  *  SP 

SO  X3(I)  =  -  X  *  SP  +  Z  *  CP 

DO  70  I=l,NMOD 
DO  65  J=1,NGRID 
X  =  FIBX(J,I) 

Z  =  FIBZ(J,I) 

FIBX(J,I)  =  X  *  CP  +  Z  *  SP 
65  FIBZ(J,I)  =-X*SP+Z*CP 

70  CONTINUE 

GO  TO  80 

3- DIHENSIONAL  COORDINATE  TRANSFORM 
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n  n  n  n  n  n  o 


C 

75  CALL  CTRN3(VB,VC,C) 

C 

DO  76  I=1,NGRID 
VO(l)  =  X1(I) 

VO(2)  =  X2{I) 

VO(3)  =  X3(I) 

CALL  TRN3{VO,VN.C) 

XI (I)  =  VN(1) 

X2(I)  =  VN’(2) 

76  X3(I)  =  Vn{3) 

DO  78  I=l,HnOD 

DO  77  J=1,UGRID 
VO(l)  =  FIBX(J,I) 

VO(2)  =  FIB'x(J.I) 

VO(3)  =  FIBZ(J,I) 

CALL  TRIJ3(VO,VN,C) 
FIBX(J,I)  =  VN(1) 
FI3Y(J,I)  =  \li{2) 

77  FIB2(J,I)  =  VN(3) 

78  CONTINUE 

WRITE  OUTPUT  FILE 

30  WRITE (21 ,110)  IJGRID,NMOD 
IF(1MAS  .E2.  0)  GO  TO  85 
DO  S3  I=1,NGRID 
63  WRITE(21,160)  I,  CMASS(I) 


85 

DO  87  I=1,NGRID 

87 

WRITE(21,210)  I,  X1(I),X2(I),X3(I) 

DO  95  I=l,NMOD 

WRITE(21 ,410)  I,  FREQ(I),  OMEGA(I) 

DO  90  J=1,NGRID 

90 

WRITE(21 ,550)  J,  FIBX(J,I),  FIBY(J,I), 
IF(ICM0D  .EQ.  0)  GO  TO  95 

DO  93  J=l,riGRID 

FIBZ(J,I) 

93 

WRITE(21,550)  J,  PSI1(J,I),  PSI2(J,I), 

PSI3(J,I) 

95 

CONTINUE 

FORMAT 

100 

FORMAT (10X,6I 5) 

104 

FORMAT(F10.4) 

106 

FORMAT(10X,3F10.4/10X,3F10.4) 

110 

FORMAT(/ 'NUMBER  OF  GRIDS  =', 15/ ' NUMBER  OF 

MODES  =' , 

150 

FORMAT(10X,2I5.D15.4) 

160 

F0R:!AT(I4,D15.6) 

200 

FORMAT(20X,3F10.3) 

210  FCRHAT(3X,I4,3X,3D15.6) 

300  FORMAT ( ax, 14) 

400  FORMAT (8X,D14. 6) 

410  FCRMAT(/ 'MODE  15/ ' FREQUENCY  = ' , FIO . 3 ,2X , ' HZ' / ' OMEGA  =',D15.6) 
500  FORMAT (2X,6D13. 6) 
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550  FORKAT(3X,I4,3X,3D15.6) 

C 

RETURN 

END 

*  MODATC:  READ  AND  TAYLOR  MODAL  DATA  OF  BODY  C  * 

SUBROUTINE  MODATC  (M,NPT ,X1 ,X2 .X3 .FIBX.FIBY ,FIBZ , 
+  CHASS, FREQ, OMEGA) 

C 

IMPLICIT  REAL*8(A-H,0-2) 

DIMENSION  FREQ(M) ,OHEGA(M) ,CMASS(NPT) 

DIMENSION  Xl(NPr) ,X2(NFT) ,X3(NPT) 

DIMENSION  FI3X(NPT,M) ,FIBy(NFT,M) ,FIB2{NPT,M) 
DIMENSION  V3(3) ,VC(3) ,VO(3) ,VN(3) ,C(3,3) 

C 

REWIND  9 
C 

PI  =  4.0  ATAN(l.DO) 

TFI  =  2.0  *  PI 
FRAD  =  Pl/iaO.O 
C 

READ (9,100)  ICOR , IMAS , ICMOD , ICTR , NGRID , NMOD 
C 

IF(ICTR  .LE.  1)  GO  TO  3 
IFdCTR  .GE.  3)  GO  TO  2 
READ (9, I 04)  ALF 
GO  TO  3 

2  READ(9,106)  V3 ( 1 ) , VB(2) , VB(3) , 

+  VC(1),VC(2),VC(3) 

C 

3  IF(IMAS  .EQ.  0)  GOTO  9 

C 

DO  5  1=1, IMAS 
READ(9,150)  J1,J2,BMAS 
DO  4  J=J1,J2 

4  CMASS(J)  =  BMAS 

5  CONTINUE 
C 

9  DO  10  1=1, NGRID 

10  READ(9,200)  XI ( I ) ,X2 (I ) ,X3(I ) 

C 

IFdCOR  .EQ.  0)  GO  TO  30 
C 

DO  20  1=1 , NGRID 
R  =  Xld) 

THTA  =  X2(I)  *  FRAD 
Xld)  =  R  *  COS  (THTA) 

20  X2(I)  =  R  *  SIN(THTA) 

C 

30  DO  50  1=1, NMOD 

READ (9, 300)  IMOD 
READ(9,400)  FREQ(IMOD) 


o  o  o  o  o  r> 


OMEGA (IMOD)  =  TPI  *  FREQ(IMOD) 

C 

DO  40  J=1,NGRID 

READ (9,500)  F I BX ( J . IMOD ) . FI BY ( J , IMOD ) , FIB2 ( J , IMOD ) 
40  CONTINUE 

50  CONTINUE 
C 

IF(ICTR  .LE.  1)  GO  TO  80 
IF(ICTR  .GE.  3)  GO  TO  75 
C 

C  2-DIHENSIONAL  COORDNATE  TRANSFORM 

C 

ALP  =  ALP  *  FRAD 
CP  =  COS  (ALP) 

SF  =  SIN(ALP) 

DO  60  I=1,NGRID 

X  =  X1(I) 

2  =  X3(I) 

X1(I)  =  X  *  CP  +  2  *  SP 
60  X3(I)  *  -  X  *  SP  +  2  *  CP 

C 

DO  70  I=l,MnOD 
DO  65  J=1,NGRID 
X  =  FI3X(J,I) 

2  =  FIB2(J,I) 

FIBX(J,I)  =  X  *  CP  +  2  *  SP 
65  FI32(J,I)  =-X*SP+2*CP 

70  CONTINUE 
GO  TO  80 

3-DlMENSIOMAL  COORDINATE  TRANSFORM 

75  CALL  CTRN3(VB,VC,C) 

DO  76  I=1,NGRID 
VO(l)  =  X1(I) 

VC(2)  =  X2(I) 

VO(3)  =  X3{I) 

CALL  TRN3(VO,VN,C) 

XI  (I)  =  Vt:(i) 

X2(I)  =  VN{2) 

76  X3(I)  *  VN(3) 

DO  73  I=l,NMOD 

DO  77  J=1,NGRID 
V0(1)  =  F1BX(J,I) 

VO(2)  =  FIBY(J,I) 

VO(3)  =  FIE2(J,I) 

CALL  TRN3(VO,VN,C) 

FIBX(J,I)  =  VN(1) 

FIBY(J,I)  =  VN(2) 

77  fIB2(J,I)  =  VN(3) 

78  CONTINUE 

WRITE  OUTPUT  FILE 
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non 


C 

80  WRITE(21,110)  NGRID.raOD 

IF(II1AS  .EQ.  0)  GO  TO  85 
DO  33  I=1,NGRID 
83  WRITE(21,160)  I,  CMASS(I) 

85  DO  87  I=1,NGRID 

87  WRITE(21,210)  I,  X1(I) ,X2(I) ,X3(I) 

C 

DO  95  I=1,MM0D 

WRITE(21,410)  I,  FREQ(I),  OMEGA(I) 

DO  90  J=1,NGRID 

90  WR:TE(21 , 550)  J,  FIBX(J,I),  FIBY(J,I),  FIBZ(J,I) 

95  CONTINUE 

FORMAT 

100  FORMAT (lOX, 61 5) 

104  FORMAT (FI 0.4) 

106  FO.RMAT(10X,3F10.4/10X,3F10.4) 

110  FORMAT(/ 'NUMBER  OF  GRIDS  ='.15/ 'NUMBER  OF  MODES  =',I5) 

150  FCRHAT(10X,2I5,D15.4) 

160  FORMAT ( 14, D1 5. 6) 

200  FORHAT(20X,F10.2,F10.4,F10.7) 

210  FORMAT(3X,I4,3:<,3D15.6) 

300  FORMAT (3X, 14) 

400  FORMAT (8X,D14. 6) 

410  FORMAT(/ 'MODE  =', 15/ ' FREQUENCY  =’ ,F10.3,2X, 'HZ' / 'OMEGA  =',D15.6) 
500  F0RMAT(24X,3D15.6) 

550  FORMAT(3X,I4,3X,3D15.6) 

r 

RETURN 

END 

*  SUBROUTINE  CTF.N3  * 

SUBROUTINE  CTRN3  (VB,VC,C) 

IMPLICIT  REAL*8(A-H.O-Z) 

INIEGER  I,I2,I3,IP2,IP3,J,K 
DIMEMSION  VB(3) ,VC(3) ,C(3,3) 

DIMENSION  IP2(3) , IP3(3) 

DI.MENSION  XM(3)  ,YN(3)  ,ZN(3) 

DATA  IP2  /2,3,1/ 

DATA  IP3  11.1.21 
C 

BL  =  0. 

DO  10  1=1,3 

10  BL  =  EL  +  VB(I)  *  VB(I) 

BL  =  DSQRT(BL) 

DO  20  1=1,3 
ZN(I)  =  VB(I)/BL 
20  C(I,3)  =  ZN(I) 

C 

DO  30  1=1,3 
12  =  IP2(I) 
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13  =  IP3(I) 

30  YN(I)  =  ZN(I2)  VC(I3)  -  2N(I3)  *  VC(I2) 

YL  =  0. 

DO  40  1=1,3 

40  YL  =  YL  +  YN(I)  *  YN(I) 

YL  =  DSQRT{YL) 

DO  50  1=1,3 
YN(I)  =  YN(I)/yL 
50  C{I,2)  =  YNd) 

C 

DO  60  1=1,3 

12  =  IP2(I) 

13  =  IP3(I) 

:;N{I)  =  YN(I2)  *  2N(I3)  -  YN{I3)  *  ZN(I2) 

60  C(I,1)  =  XN(I) 

RE  TORI'! 

END 

*  SUBROUTINE  TRI'IS  * 

A  A  A  A  A  A  A  AA  A  A  AAA  A  AAAAA  AAAA  A  A  A  A  A  A  A  A  AA 

SUBROUTINE  TR!'13  (VO,Vn,C) 

IMPLICIT  REAL*8(A-H,0-Z) 
ir.'TEGER  I,J,K 

DIMENSION  VO(3),VN(3),C(3,3) 

C 

DO  10  1=1,3 
VN(I)  =  0. 

DO  10  J=l,3 

10  VN(I)  =  VN(I)  +  C(1,J)  *  VO(J) 

RETURl! 

END 
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C  CALCULATE  A  MATRIX  (  UDOT  COEFFICIENTS) 

FIXED  lER, IPVT,NPT,M,N,NPTC,HC 
PARAH  MPT=17  ,MPTC  =^63,H=4,MC=4 
D  DIMENSION  A(9,9) 

D  DIMENSION  U1D21 ( 10 , 10 ) , U1D22( 10 , 10 ) , U1D32 ( 10 , 10 ) 

D  DIMENSION  UJD21 ( 10 , 10 ) , UKD21 ( 10 , 10 ) , BU222 ( 10 , 10 ) , BU223 ( 10 . 10 ) 

D  DIMENSION  B22UUJ ( 10 , 10) , B22UUK( 10 , 10) , B32UUJ ( 10 , 10 ) , U1D31 { 10 , 10 ) 
D  DIMENSION  B3U213 (10 , 10 ) , B32U»'K(10 , 10 ) 

D  DIMENSION  FIBX( 17 , 10 ) , FIBY{ 1 / , 10 ) , FIBZ ( 17 , 10 ) , PSI 1 { 10 ) 

D  DIMENSION  FICX(63 , 10 ) , FICY( 63 , 10 ) , FICZ (63 , 10 ) , PSI2 ( 10 ) 

D  DIMENSION  PSI3 ( 10 ) , XI (63 ) , X2(63 ) , X3 (63 ) 

ARRAY  IFVT(20) 

ARRAY  U1D12(10) , BU221 ( 10 ) ,U1D13 ( 10 ) ,UJD11(10) ,UKD11(10) 

ARRAY  B12UUJ(10) , B12UUK(10) ,B3U211(10) 

ARRAY  STRNB(IO) ,STRNC(10) ,OMGAB(10) ,OMGAC(10) 

INTEGER  I,J,K,L 

FARAM  C1=.4,C2=10. ,C3=24.65,DELT=.010 
INITIAL 

*  READ  IN  VALUES  OF  CONSTANTS 
READ  (25,107)IDUM 
READ  (25,102)  UlDll 
READ  (25,107)IDUM 
DO  1  J=  1,10 

READ(25,I03)L,U1D12(J) ,UJD11(J) ,B12UUJ(J) ,BU221(J) ,OHGAB(J) 

1  CONTINUE 

READ  (25,107)IDUM 
DO  2  K  =  1,10 

READ(25,103)L,B12UUK(K) ,U1D13 (K) , UKDll (K) ,B3U211(K) ,OMGAC(K) 

2  CONTINUE 

READ  (25,107)IDUH 
DO  3  I  =  1,10 

DO  4  J  =  1,10 

READ(25,106)L,L,UJD21(J,I) ,UiD21(J,I) ,BU222(J,I) ,B22UUJ(J,I) 

4  CONTINUE 

3  CONTINUE 

READ  {25,107)IDUH 
DO  5  K  =  1,10 
DO  6  I  =  1  , 10 

READ(25,104)L,L,U1D32(I,K) ,B3U213(I,K) ,B32UUK(I,K) 

6  CO.MTINUE 

5  CONTINUE 

READ  (25,107)IDUM 
DO  7  K  =1,10 
DO  3  J  =  1,10 

READ(25,106)L,L,U1D22(J,K) ,UKD21(J,K) ,BU223(J,K) ,B22UUK(J,K) 

3  CONTINUE 

'!  CONTINUE 

READ  (25,107)IDUM 
DO  98  K  =  1,10 
DO  99  J  =  1,10 

READ  (25,105)  L , L , U1D31 ( J , K) , B32UUJ ( J , K ) 

99  CONTI  in 'Z 

93  CONTINUE 

1C2  FORMAT  (D15.6) 
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103  F'^RtlAT  (I4,5D13.6) 

104  FCRMAT(2I4,3D15.6) 

105  FORMAT  (214,2015.6) 

106  FORMAT  (214,4015.6) 

107  FORMAT  (15) 

DO  94  J  =  1,10 
DO  95  I  =1,MPT 

READ(25,900)FIBX(I,J),FIBY(I, J),FIBZ(I,J) 

95  COMTIMUE 

94  COMTIMUE 

DO  96  J  =  1,10 
DO  97  I  =  1,N?TC 

RE.AD  (25,900)  FICX(  I ,  J )  ,  FICY(  I ,  J )  ,  FIC2(  I ,  J ) 
97  CONTINUE 

96  CONTINUE 

DO  93  J  =  1,10 

READ  (25,900)  FSI 1 ( J ) , PSI2 ( J ) , PSI3 ( J ) 

93  CONTINUE 

DO  92  I  =1,NFTC 

READ  (25,900)  XI ( I ) , X2 ( I ) ,X3 ( I ) 

92  CONTINUE 

900  FORMAT  (3D20.6) 

CALL  STRNEG  (M , CMGAB , STRNB ) 

CALL  STRNEG  (HC , OMGAC , STRNC) 

DO  91  J  =1,H 

W.RITE  (21,900)  FSIl(J)  ,PSI2(J).PSI3(J) 

91  CONTINUE 

DERIVATIVE 

NOSORT 

TA  ^  10.-10,*STEP(6.35) 

T1  =  2. 

TAOl  =  80.*TIHE**2 
TA02  =  320.*(1.-STEP(5.  )) 

TA  =  3WITCK(TIME  .LE.  T1 , TAOl . TA02 ) 

*  TA  =  5. 

TA  +  10.-10.'^STEP(6.85) 

AllQB  =  0. 

B12UJ  ^  0. 

DO  40  J  =  1  ,H 

AllQB  =  Q(J+1)  *  U1D12(J)  +  AllQB 

B12UJ  =  B12UUJ(J)  -  U(l+J)+  B12UJ 
DO  10  I  =J,M 

A(1*J,1  +  I)  =  UJD21( J, I  ) 

A(l  +  I , 1+J)  =  A(l+J,  l  +  I  ) 

10  CCNTINUE 

AIJIQB  =  0. 

BJU2QB  =  0. 

BJ2UUJ  =  0. 

DO  DO  I  =  1,M 

AlJlOB  =U1D21(J,I)’Q'1+:)  +  AIJIQB 
DJU2QB=  BU222(  J  ,  I ) 'Qa  +  : )  5JU2QB 
BJ2UUJ=B22UUJ( J, I l^Uf I+l ;  +  BJ2UUJ 
CONTINUE 
AlJiOC  =  0. 
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BJU2QC  =  0. 

BJ2UL'K  =  0. 

DO  30  K  =  1,I!C 

A1J10C=  U1D22(J,K)  *Q(l+tI+K)  +  AIJIQC 
A(1+I1+K,1  +  J)=UKD21(J,K) 

A(l  +  J,  1+M-"K)  =  Ad+M+K,  1+J) 

BJU2QC=  BU223(J,K)'^Q;i+tt+K)  +  BJU2QC 
BJ2UUK  =  B22UUK(J,K)  *U(l+M-t-K) 

30  COtJTItIUE 

A(l+J,l)  =  UJDll(J)  +  AIJIQB  +  AIJIQC 
A(l,l+J)  =  A(l+J,l) 

B(1+J)=U(1)*U(1)^(B'J221(J)+BJ'J2QB+BJU2QC)  .  .  . 

+2*'U(1)*(BJ2UUJ  +  BJ2UUK)-  STRMB ( J ) *Q ( 1+ J ) 

40  CONTINUE 

AllQC  =  0. 

312UK  =  0. 

DO  60  K  =  1,HC 

B12UK  =  B12UUK(K)*U(1+M+K)  +  B12UK 
AIKIQC  =  0. 

AIKIQB  =  0. 

BKU2QC  =  0. 

BK2UUK  =  0. 

DO  50  I  =  1,MC 

AIKIQC  =  0(1+M+I)  *  U1D32(I.K)  +  AIKIQC 
EKU2QC  =  B3U213(I,K)*Q(1+N+I)  +  BKU2QC 
BK2UUK  =B32UUKa.K)*Ud+M-I)  +  BK2UUK 
A( l+H+I , i+n+K)  =  0. 

50  CONTINUE 

PKU2QB  =  0. 

BK2UUJ  =  0. 

DO  5  5  J  =  1  ,M 

AIKIQB  =  Oa+Jd  U1D31(J,K'  +  AIKIQB 
3KU203  =  BU223(J,K)  *Q;i+J)  +  EKU2QB 
DK2UUJ  =B32UUJ(J,K)*Ui 1+J;  +  BK2UUJ 
5  5  COt-'TINUE 

AllQC  =  QC+N+K;  *  U1D13(K'  +  AllQC 
A^  l  +  il+K,  1+M+K)  =  1 . 

A(1+!H-K,1)  =  UKDli(K)  +  AIKIQB  +  AIKIQC 
Ad,  1+I1+K)  =  Ad+N+K,l) 

B'i*lI-^K)  =  Ud)*Ud)^(B?-J2il(K)  +  BKU2QB  +  BKU2QC)+  ... 

D’-U  (1 )  *  ( BK2UUJ  +  L”2'J-JK )  -STRNC  { K  )  Q  { 1+tl+K  ) 

50  'TCNTINUE 

Ad,i)=  UlDil  +  AllQB  +  AllQC 

5  1}  =  2'Ud)*(  B1..UJ  +  B12UK  '•+  TA 

'.•DIT£(22  .ZOOM  (A(  I  ,  j;  ,  J^l  ,1*:UHC)  .  1  =  1  ,  lt-11+MC) 

ZOO  r:?.:iArd‘Ei4.6) 


ALI.  ZCEFA  A  ,  9 . 9  ,  IF7T  .  lEF,' 


•■:?ITE(8 , 300  )  lE.R 
r^PNAT  dZO) 

tf'IEF.::e.  o)  go  to  loo 

CALL  DGESL  ( A , 9 , 9 . IFVT . B , 0  ) 


v:rite'4,2''‘C;  dd' 
■J  =ii;tgpz'  0. 


■  =1 ,  i+nmc) 
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Q  =  INTGRL(0. ,U,9) 

'CALCULATE  VARIOUS  DISPLACEMENTS  AT  THREE  LOCATIONS  AND  ROTATIONS  AT  TIP 

:-:hdi3p  =  o. 

VKDISP  =  0. 

ZHDISF  =  C. 
aHROT  =  0. 

YHROT  =  0. 

ZHROr  =  0. 

DO  70  I  =  1,M 


DO  70  I  =  1,M 

aHDISP  =  :<HDISP 
YKDISP  =  YHDISP 
ZHDISP  =  Zr'DISP 
XHROr  =  XHROT  + 
YHROT  =  YHROT  + 
2HR0T  =  2HR0T  + 
CONTINUE 
aIREL  =  0. 

YIREL  =  0. 

ZIREL  =  C. 

Z4REL  =  0. 

Y4FEL  =  0. 

24REL  =  0. 

DO  80  I  =  1,HC 

aIREL  =  XIREL  + 
YIREL  =  YIREL 
ZIREL  =  ZIREL  + 
a4REL  =  :':4REL  + 
Y4REL  =  Y4REL  + 
Z4REL  =  Z4REL  + 
COHTIIIUE 

AlCISP  =  aHDISP  +  5 


+  Q(I+1)*FIBX(NPT,I) 
+  Q(I+1)'FIBY(NPT,I) 
+  Q(I+1)*FIBZ(NPT,I) 
Q(I  +  1)’'PSI1(I) 

Qd  +  D'PSIZd) 
Q(I*1)*PSI3(I) 


Q{i+i+ri;*ric:<(i,i) 

Q(I  +  1+M)-‘FICV(1 , 1 ) 
Q(I+1+M)*FICZ(1 , I ) 
Q(I  +  l+MdFICX(4, 1) 
Q(I+1+M)-FICY(4,I) 
q;I+1+M)*FIC2(4,I) 


aHDISP 

VICISP  =  YHDISP  +  YIREL  +  (ZHROT^Kl(l)  -  XHR0T*A3(i)) 

ZIDISP  =  ZHDISP  +  ZIREL  +  'XHROT*X2ll)  *  VHROT^XKU) 

a4PIS?  =  XHDISP  +  X4FEL  +  fYHROT*X3(4)  -  ZHROT^X2(4;; 

YiDISP  =  YHDISP  +  Y4REL  +  ( ZHROT ''Xl  ( 4 )  -  XrROT*X3(4)) 

Z;DIS?  =  ZHDISP  +  Z4REL  -*•  (XHROT’X2(4)  -  YHROr'Xn  4 )  ) 

RETURN 

100  WRITE',  6,101)  TIME,  lER 

101  F'F:IAT(  '0  lER  =M7) 

TAIL  EI!DJOB 

PRINT  TA,5(1-?) ,U(l-5) 

C'H.'TRL  F II;TIM=10.  .DELPRT  =  .1 

SAVE  . 0 1 5 , TA , Q ( 1  )  , 0 ( 2 ) , Q ( 3 ) , L ( 4 1 . Q ; 5 ) . Q { 6 ) , 5 ( 7 ) , A I  3  d  ? ( 9 ) _ 

'V'  1  )  ,XKDISF,  YHDISP  ,  ZHDISP,  XHF/^r,  YHROT,  ZHROT  ,X1D  I  SR  .... 

Y 1 D I SP , Z 1 DI S P , X4DI 3P , Y4D I S P , Z4r  I SF 
'■PAIH  '  01  ,:il  =  7  ,DE=TEK618;  TIXECII^F  .LE=10  .  ,UM=  ... 

SEC),  TA' UN=  '  LB-IN'  ) 

'?A?H  'OZ.NI^?  ,DE  =  TEK618)  TINE  '  N I  =  5  .  LE- iO  .  ,  U:r=  '  SEC  '  ... 

"  ,v  1 ' 'u:;='PAD' / 

‘FAFH  '03,;il  =  7  ,DE^TEKG15)  TINE  aiI-5  .  LE^IO  .,  UN'^  '  SEC  ;  ... 

.V'.  i:(ui;=  rad;  SEC; 

'P.AFH  '  :i,NI  =  7  ,DE=rEK618)  TINE  (NI  =  5  .  LE^IO  ..  UN=  '  SEC  .)  ... 

.'/C  IN'  ; 

•P.AFH  'G5,:;i  =  7,DE  =  rEK6ie)  TINE'NI^F  .LE^lj  .  , 'JN=  '  .SEC  )  ... 


XIREL 

YIREL 

ZIREL 

X4FEL 

Y4REL 

Z4REL 


{YHRCT*X3(1) 
(ZHROT^Xl(l) 
';XHROT*X2ll ) 
fYHROT*X3(4) 
(ZHROT "Xl (4) 
(XHROT’X2(4) 


ZHROT^X2( 1 ) ) 
XHR0T*A3(i)) 
YHROT* XK  1  ^  ) 
ZHROT* X2 (4; ; 
XrROT*X3(4) ) 
YHR0r*AC4  )  ) 


‘FAFH 


TINE(NI=5 , 


UN= ' 3E' 


UN^'SEC 


*  ,  Q(3) (UM=' IM' ) 

*RAFH  (G6 .MI=7 ,DE^TEKG18)  TIME (MI-8 , LE^ 10 UM= ' SEC ' )  ... 

*  ,Q(4) ) 

^RAPH  (G7,MI=7,DE=TEK618)  TII1E(NI=5 , LE-10 , , UM= ' SEC ' )  ... 

*  ,Q(5)(UM='1N' ) 

*RAFH  (G7A,NI-7,DE=TEK618)  TIME (MI=5 , LE=10 ., UN= ' SEC ' )  .. 

*  ,Q(6) (UN=' IM' ) 

■^RAFH  (G7B,MI=7  ,DE=TEK613)  TIME  (NI-5  ,  LE=10  .,  UN=  '  SEC  '  )  .. 
^  ,Q(7) (UN=‘ IN' ) 

-r.^FH  (G7C,MI  =  7  .DE-TEK613)  TIME  (NI  =  5  ,  LE-10  .,  UH=  '  SEC  '  )  .. 

*  ,Q(3)(UN='IM' ) 

*RAPH  (G7D,MI=7,DE=TEK613)  TIKE (NI=5 , LE=1 C ., UN= ' SEC ' )  .. 

*  ,  Q(9) (UM=' IN' ) 

*RAFH  (G8,MI=7 ,DE=TEK&18)  TIME (MI=5 , LE-10 ., UN= ' SEC ' )  ... 
^  .:-:HDISP(UM=' IN' , LI-5, LO—19.  ,SC=2.  ) 

GRAP.M  (G9,NI  =  7  .DE-TEK618}  TIME  (NI-5  ,  LE-10  .  ,  UN- '  SEC  '  )  ... 

,YHDISF(UM='  IN'  ,  LI-4,  L0--2.  .SC-1 . 

*.RAFH  (G.A,ni-7  ,DE-TEK618)  TIME  (i;i-5  ,  LE-10  .,  UN- '  SEC  '  )  ... 
,ZHDI5P{UN-' IN' , LI-4) 

GRAPH  (GR3,MI-7  ,DE=TEK618)  TIME (MI-5 . LE-10 ., UN- ' SEC ' )  .. 

, MHROT ( UN- ' RAD ' , LI-4 . L0=- . 020 , SC- . 005 ) 

^.RAPH  (GP,9  ,MI  =  7  ,DE=TEK618)  TIME ( MI -5  , LE-i 0  .  ,  UN- '  SEC  '  )  .. 

*  ,VHROT(U:i=  '  PAD'  ) 

'RAPH  (GRA,NI=7 ,DE-TEK618)  TIME (MI- 5 , LE-10 . , UN- ' SEC ' )  .. 
,  .CHRCT  '  UM-  '  RAD  '  ) 

’RAFH  (GB, MI-7, DE=^TEK81S)  TIME  (MI-5  ,  LE- 1 0  .,  UM- '  SEC  '  )  ... 

*  .••:iBISFfUn='  lil'  ,LI  =  5.L0— 2?.  ,SC  =  5  .  ( 

G.c.ARH  'GC,M:-7  ,rE  =  TEK61S)  TIME  <  r;i  =  5  ,  LE-10  .,  UM- '  SEC  '  )  ... 

,ViD;SF(UM=' IM' ,LI=4.L:=-3. ,SC=2. ) 

’GRAPH  (GD,;iI=7 ,DE=TEK6i8)  TIME ; MI-5 , LE-10 ., UM- ’ SEC ' )  .. 
ziDiSF'UM- ' i:r  ) 

•GRAFR;  (  GE  ,MI-:  .  DE-TEK61S)  TinE-NI-5.LE  =  l?.  ,UM=  SEC  )  .. 
"  ,:-:4DisPsU:;=' IN' ,Li=  5  ,lo--3.  s:-:.  > 

GRAPH  'GF,MI-7 .DE-TEK613'  TIME : Ml =5 , LE-1 0 . , UN- ' SEC ' )  ... 

,  Y4DISF  ;UM=  '  IN'  .  LI-4  ,LC  =  -1  .  I ,  S'.'-l  .  , 

"GRAPH  ('GG,';!-' .DE-TEKGie)  PINE  '  MI  -  5  ,  LE-i  0  .  .  UM- '  SEC  '  '  .. 

,C4DIC'F  'CM-' IN'  ) 

LABEL  fGl)  APPLIED  TOF-CUE 

l'BEl  f  :-2 !  a;;gular  displacement 

LLSEL  (Gi)  ANGULAR  VELOCITY 
L.AB-L  '  4;  -EriERALIZED  DISPLA'CNENT  Cl 

LABEL  'GSi  GEMEFALIZED  DISPLA'IEMCNT  .’2 

LAPEL  n;G;  GENERALIZED  DISPLA:FNEN:  ,3 
LAPEL  ' G7  )  GENERALIZED  D I SF LA" Eirii I  .4 
LAPEL  ,  G7A)  GEM’ERALIZED  CIGF  LA  .'FN-Nr  C5 

LA-PL  Gl’H;  GEMEFALIZED  DI5FLA  EN'NT  .A 

lAGFL  'G7C)  GEMEFALIZED  LCrLA'PNlNT 

'APFL  '  '  G-NFRALIZED  IIS^LCFNfNC  ,I 

LAPEL  33,  lEFLECTI  M  AC  N  T'E 

1  ■  1  -  "  -  • " ‘  1  :  E  L  E "  T I ' '  N  AT  BN  N  Tl.- 

.'A  -  Z  i  “rLECTI  .M  AT  F  !'  .  .  .* 

1A2EL  ’  ■?-)  FNIATI  N  AT  ?'  N  TIE 

lAPEL  '  T '-*  )  F'.TATN'M  AT  P  '  N  LIE 
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LABEL  (GB)  X  DEFLECTION  AT  DISH  POINT  1 

LABEL  (GC)  Y  DEFLECTION  AT  DISH  POINT  1 

LABEL  (GD)  Z  DEFLECTION  AT  DISH  POINT  1 

LABEL  (GE)  X  DEFLECTION  AT  DISH  POINT  4 

LABEL  (GF)  Y  DEFLECTION  AT  DISH  POINT  4 

LABEL  (GG)  Z  DEFLECTION  AT  DISH  POINT  4 

END 
STOP 
FORTRAN 

*  STRAIN  Et:ERGY 

SUBROUTINE  STRNEG  (M, OMEGA, STRAIN) 
REALMS  OMEGA(IO) ,STRAIN(10) 

DC  10  I  =  i,n 

STRAIM(I)  =  OnEGA(I)**2 
10  CONTINUE 

RETURN 
END 


APPENDIX  F 


FIGURES 
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Figure  A.  1 


Undeformed  Reflector 


Figure  A.2 

First  Mode  of  Reflector 
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Figure  A.3 

Second  Mode  of  Reflector 
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Figure  A. 5 

Fourth  Mode  of  Reflector 
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MICROCOPY  RESOLUTION  TEST  CHART 

NATK)NAl  BUREAU  OF  STANDARDS -1963- A 


0% 


S  X  5  X  ^ 
•  •  •  •  •  • 
in  ^  ^  ^ 


in 


T^xTnrarrnn 


Figure  B.  1 

Angular  Velocity  of  Stiffer  Reflector  at  -155* 
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Tire  tSECI 


Figure  B.3 

Fourth  Boom  Generalized  Coordinate  of  Stlffer  Reflector  at  -155* 
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Figure  B.5 

Third  Reflector  Generalized  Coordinate  of 
Stiffer  Reflector  at  -155* 
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Figure  B.6 

Third  Reflector  Generalized  Coordinate  of 
More  Flexible  Reflector  at  >155* 
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Figure  B.8 

Vertical  Deflection  at  Boom  Tip  of 
More  Flexible  Reflector  at  -155* 
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Figure  B.9 

Rotation  of  Boom  Tip  About  Vertical  Axis 
of  Stlffer  Reflector  at  -153* 
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Figure  B.IO 

Rotation  of  Boom  Tip  About  Vertical  Axis 
of  More  Flexible  Reflector  at  -155* 
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Figure  C.l 

Angular  Displacement  of  Stiffer  Reflector  at  -135* 
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Figure  C.3 

Angular  Displacement  of  Stiffer  Reflector  at  -155* 
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Figure  C.4 

Horizontal  Deflection  of  Boom  Tip  of  Stiffer  Reflector  at  -135' 
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Figure  C.5 

Horizontal  Deflection  of  Boom  Tip  of  Stiffer  Reflector  at  -145* 
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Figure  C.7 

Rotation  of  Boom  Tip  About  Horizontal  Axis 
of  Stiffer  Reflector  at  >135' 
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Figure  C.8 

Rotation  of  Boom  Tip  About  Hoiisontal  Axis 
of  Stiffer  Reflector  at  -145* 
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Figure  C.9 

Rotation  of  Boom  Tip  About  Horizontal  Axis 
of  Stiffer  Reflector  at  -155* 
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Figure  C.12 

Horizontal  Deflection  of  DUh  Point  1  of  Stiffer  Reflector  at  -135* 
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Figure  D.  1 

First  Boom  Generalized  Coordinate  of  Stiffer  Reflector  at  >155* 
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Figure  D.2 

First  Boom  Generalized  Coordinate  of 
Stiffer  Reflector  Tilted  5*  Out  of  Plane 
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Figure  D.3 

Fourth  Reflector  Generalized  Coordinate 
of  suffer  Reflector  at  -155* 
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Figure  D.4 

Fourth  Reflector  Generalized  Coordinate 
of  Stlffer  Reflector  Tilted  5*  Out  of  Plane 
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Figure  D.5 

Out-of-Plane  Deflection  at  Boom  Tip  of  Stiffer  Reflector  at  -155* 
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Figure  D.6 

Out-of'Plane  Deflection  at  Boom  Tip  of 
Stiffer  Reflector  Tilted  5*  Out  of  Plane 
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Figure  D.7 

Out-of-Plane  Deflection  at  Dish  Point  4 
of  Stiffer  Reflector  at  -155* 
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Figure  D.8 

Out-of'Plane  Deflection  at  Dish  Point  4 
of  suffer  Reflector  Tilted  5*  Out  of  Plane 
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